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Einführung 


BASIC ist heute so etwas wie eine Standard-Programmiersprache, und für viele private 
und berufliche Computer-Verwender ist BASIC die als erste gelernte Sprache. Das liegt 
vor allem an der leichten Erlernbarkeit — selbst Anfänger können bereits nach wenigen 
Stunden eigene Programme schreiben. Praktische Gründe für die Dominanz von BASIC 
sind aber auch die Dialogfähigkeit, die interaktives Arbeiten am Computer ermöglicht 
(Mensch-Maschine-Dialog), und die Tatsache, daß die am meisten verbreiteten Tischcom- 
puter nur BASIC „verstehen“. 

Häufig hört man von erfahrenen Programmierern, BASIC weise eine Reihe ernst zu 
nehmender Mängel und Nachteile auf, vor allem 


— BASIC erlaubt keine strukturierte Programmierung, verleitet vielmehr zu extensiven 
Verzweigungen und Verschachtelungen (,,Spaghetti’-Stil); 

— es sind nur globale Variablen möglich (nur solche, die für das ganze Programm gelten), 
lokale Variablen z. B. in Unterprogrammen können nicht definiert werden; 

— symbolische Adressierung ist in der Regel nicht möglich, Variablennamen sind oft viel 
zu kurz (meist nur 2 Zeichen); 

— das Aneinanderhängen mehrerer Programme (chain, append oder merge) durch Zu- 
laden vom Massenspeicher ist normalerweise nicht möglich; 

— WHILE-Schleifen, CASE-Strukturen und IF-THEN-ELSE werden nur äußerst selten 
geboten. 


Das stimmt natürlich alles. Es sollte trotzdem immer ehrlich abgewogen werden, ob man 
auf all diese fehlenden Möglichkeiten nicht auch verzichten kann. Die weite Verbreitung 
von BASIC läßt den Schluß zu, daß in sehr vielen Anwendungsfällen tatsächlich darauf 
verzichtet werden kann. 

Es gibt jedoch in professionellen Anwendungsbereichen sehr schnell Fälle, die den Pro- 
grammierer sofort mit den aufgeführten Mängeln konfrontieren. Das wissen natürlich 
auch Computerhersteller. Und so gibt es inzwischen einige Tischcomputer auch unter 
10 000,— Mark, die mit erweiterten BASIC-Versionen arbeiten, oft als enhanced oder 
extended BASIC bezeichnet. Manchmal werden sogar BASIC-Dialekte angeboten, die 
eine Quasi-Strukturierung zumindest beim Ausdrucken der Anweisungslisten ermög- 
lichen. 

Auf der anderen Seite wird man zunehmend mit Pascal konfrontiert. Und es kann gesagt 
werden, daß die Programmiersprache Pascal die eben bei BASIC aufgezeigten Mängel 
nicht aufweist. Es gibt aber andere Nachteile: 


— Pascal ist nicht ganz so leicht erlernbar wie BASIC; 

— interaktives Arbeiten ist nicht so ohne weiteres möglich, weil Pascal kompilierend 
umgesetzt wird; 

— die äußerst bequemen Möglichkeiten von Stringmanipulationen in der Programmier- 
sprache BASIC existieren in Pascal nicht; 

— Ein-/Ausgabeoperationen werden von der Sprachendefinition her nicht unterstützt. 


Einführung 


Es gibt aber entscheidende Vorteile in Pascal. Vor allem zu nennen sind die Möglichkeiten, 
verschiedene Variablentypen mit relativ „langen’’ Variablennamen zu definieren, und die 
Tatsache, daß man sozusagen zu strukturierter Programmierung gezwungen wird. 


In diesem Band der Programmbibliothek haben wir eine Reihe von Beiträgen zusammen- 
gestellt, in denen BASIC und Pascal in sehr eingängiger Weise gegenübergestellt sind. Vor- 
angestellt ist der Aufsatz „Graphische Darstellung von Programmabläufen‘’ von Wolfgang 
Schneider. Darin werden die beiden Möglichkeiten der Programmbeschreibung bzw. 
-dokumentierung diskutiert, die jeweils gerade den beiden hier verglichenen Programmier- 
sprachen angemessen erscheinen: die konventionellen Programmablaufpläne für „ältere’' 
Sprachen und Struktogramme für Pascal. 

Es folgen zwei Beiträge von Rüdeger Baumann mit algorithmisch betrachteten Spielen, 
worin die „Vorzüge der Programmiersprache Pascal gegenüber BASIC’’ besonders gut 
herausgestellt werden können. Es sei an dieser Stelle der durch Rüdeger Baumann dem 
Herausgeber übermittelte Stoßseufzer von Sidonie Trampler (Literaturzitat [17, S. 104] 
im Beitrag „Ein Spiel — ...’‘) angeführt: „Der Herr behüte uns vor den von Didaktikern 
zum Zweck des Mathematiklernens erfundenen ‚Spielen’!’‘. Die Spiele für den BASIC/ 
Pascal-Vergleich sind gerade nicht von dieser Sorte. 

Ein großer Block mit 10 Beiträgen stammt von Dietmar Herrmann. Die einzelnen Bei- 
spiele werden anschaulich eingeführt, und es wird gezeigt, daß die verwendeten Algorith- 
men meist direkt in Pascal programmiert werden können. Für Lösungen mit BASIC müs- 
sen häufig umständliche Simulierungen und Berechnungsroutinen geschrieben werden, 
die sinnvollerweise in der Regel als Unterprogramme aufzurufen sind. 

Der letzte Block in diesem Bibliotheksband stammt von Karl Achilles. Anhand von fünf 
Beispielen wird auch dabei aufgezeigt, wo jeweils Vor- oder Nachteile bei der Programmie- 
rung in BASIC bzw. Pascal zu sehen sind. Auch hierbei wird deutlich, daß die aufgestell- 
ten Algorithmen zumeist direkt eine entsprechende Pascal-Struktur widerspiegeln. 
Möglicherweise wird durch solche Veröffentlichungen, wie die hier vorliegende, der 
Wunsch nach Pascal verstärkt. Es sollte aber nicht vergessen werden, daß BASIC bei An- 
fängern leichter die Schwellenangst abbauen kann bzw., andersherum, daß durch erst- 
malige Konfrontation mit Pascal manchmal Abneigungen verstärkt werden können. Für 
alle in der Prozeßdatenverarbeitung tätigen Anwender bleibt Pascal wegen der fehlenden 
Ein-/Ausgabe-Strukturen sowieso zweite Wahl, es sei denn, die Assembler-Sprache wird 
zur Programmierung von Schnittstellentreibern akzeptiert. Wahrscheinlich aber wird 
dieser spezielle Anwenderkreis eher zu „‚prozeßgeeigneten’’ Sprachen tendieren und dann 
auf beispielsweise Forth stoßen, eine wenigstens 10 Jahre „alte’’ Sprache, die aber erst in 
jüngster Zeit auch für Mikrocomputer interessant wird. Oder man gelangt zur neuen 
„‚Universalsprache‘’ Ada, die vom US-Verteidigungsministerium (Department of Defense, 
DoD) kontrolliert und lizenziert wird. (Ada-Compiler dürfen nur so heißen, wenn die 
DoD-Lizenz vorliegt!) 

Wir bleiben aber zunächst noch in der praktischen Gegenwart und stellen im folgenden 
BASIC und Pascal im Vergleich gegenüber. 
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Graphische Darstellung von Programmabläufen 


von Wolfgang Schneider 


Vor der Programmierung eines Problems in einer beliebigen Pro- 


grammiersprache empfiehlt es sich, 


- das Problem aufzubereiten und 


- den Programmablauf graphisch darzustellen. 


Erst anschließend sollte man, zumindest bei umfangreichen Pro- 
blemen, zum Schreiben des Primärprogrammes in einer beliebi- 


gen Programmiersprache übergehen. 


1 _PROBLEMAUFBEREITUNG 
Zur Problemaufbereitung gehört 


- eine vollständige Formulierung der Aufgabe und 


- eine Problemanalyse der Aufgabe. 


Die Aufgabe ist zunächst vollständig mit allen Randbedingungen 
in der Umgangssprache zu formulieren. In der darauf folgenden 


Problemanalyse ist u.a. zu untersuchen, 


- ob die Aufgabe überhaupt mit Hilfe einer Datenverarbeitungs- 
anlage (kurz DVA) gelöst werden kann, 

- welche alternativen Lösungswege sich für die Aufgabe anbie- 
ten und 


- welcher der möglichen Lösungswege der günstigste ist. 


2 DARSTELLUNG VON PROGRAMMABLÄUFEN 


Nachdem bei der Problemaufbereitung ein günstig erscheinender 
Lösungsweg gefunden wurde, empfiehlt es sich vielfach, die 


einzelnen Schritte zur Lösung des Problems graphisch darzustel- 
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len. Hier stehen dem Programmierer zwei wichtige Darstellungs- 
weisen zur Verfügung: 
— Darstellung mit Hilfe von Programmablaufplänen. 


— Darstellung mit Hilfe von Struktogrammen. 


Diese beiden Darstellungsweisen sollen im folgenden besprochen 


und miteinander verglichen werden. 


allgemeine schließe 
Operation Ventil A 


Verzweigung 


Ventil A 
geschlossen? 


Bemerkung: Der Text muß eine Frage (Bedingung) enthalten, die entweder 
mit ja oder nein zu beantworten ist. Je nach Beantwortung der Frage wird 
das Programm mit dem „‚Ja’- oder „Nein’’-Zweig fortgesetzt. 


Ausgabe 


Grenzstelle 


STOP START 


Bemerkung: Die Grenzstelle ist das Sinnbild für den Beginn oder das Ende 
eines Programmes 


Übergangsstelle 


Bemerkung: Falls ein Programmablaufplan auf einem anderen Blatt fortge- 
setzt werden muß, kennzeichnen die Übergangsstellen mit gleichen Zahlen, 
an welcher Stelle das Programm auf dem anderen Blatt fortgesetzt werden 
muß. 


Richtungspfeil Bemerkung: Falls keine Richtungspfeile angegeben wer- 
den, wird der Programmablaufplan von oben nach unten 
und von links nach rechts gelesen. 





Bi ld 1: Nach DIN 66001 genormte Sinnbilder von Programm- 


ablaufplänen: 


4 


Graphische Darstellung von Programmabläufen 


2.1 Programmablaufpläne 


Ein Programmablaufplan ist eine graphische Darstellung, die 
den Arbeitsablauf einer Problemstellung in einzelnen Schritten 
darstellt. 


Die Programmablaufpläne setzen sich aus verschiedenen Sinnbil- 
dern zusammen, die nach DIN 66001 genormt sind. Durch Einfügen 
eines Textes in die Sinnbilder wird die Art der Vorgänge genau 
spezifiziert. Die Reihenfolge der Vorgänge wird durch Pfeile 
angedeutet. Für einfache Aufgaben genügt die Kenntnis der in 
Bi 1 d 1 aufgeführten Sinnbilder. 


Das Zeichnen der Programmablaufpläne wird durch Schablonen 
erleichtert (B i1d 2). 





Bi ld 2: Beispiel einer Schablone zur Erstellung von Pro- 


grammablaufplänen 


Vorteile bei der Anwendung von Programmablaufplänen: 


Der Programmablaufplan erweist sich bei umfangreichen Aufgaben 
als sehr zweckmäßig. Insbesondere sind folgende Vorteile zu 


nennen: 


o Der Programmablaufplan verschafft dem Programmierer durch die 
logische Gliederung einen Überblick über den Gang der Rech- 


nung. 
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o Der Programmablaufplan verhindert das Programmieren von 
"Sackgassen". 

Eine Sackgasse würde sich z.B. in einem Programm ergeben, 
wenn ein Zweig einer Verzweigung durch Vergeßlichkeit des 
Programmierers nicht weiter berücksichtigt würde. Wenn bei 
einer späteren Benutzung des Programms dieser Zweig gewählt 
wird, so gibt es keine darauffolgende Anweisung. Das Pro- 
gramm endet in einer Sackgasse. In einem Programmablaufplan 
sind derartige Sackgassen gut zu erkennen und können somit 
vermieden werden. 

o Der Programmablaufplan bietet ein gutes Verständigungsmittel 
zwischen einem Spezialisten und einem Programmierer. 
Spezialisten verfügen teilweise über keine ausreichenden 
Programmierkenntnisse. Ein Programmierer hingegen verfügt 
nicht immer über die notwendigen Spezialkenntnisse, um pro- 
grammierbare Regeln aus einer Aufgabenstellung abzuleiten. 
Hier bietet 3ich der Programmablaufplan als gemeinsames 
Verständigungsmittel an. 

o Der Programmablaufplan hilft bei der Fehlersuche von logi- 


schen Fehlern. 


Durch dıe logische Gliederung des Problems in eine Folge von 
Einzelschritten ist der Programmablaufplan wegen der besse- 

ren Übersicht meist besser als das Programm selbst geeignet, 
logische Fehler im Programmablauf zu finden. 


Der Programmablaufplan dient zur Dokumentation des Programms. 


Programme sollen auch später, eventuell von anderen Personen, 
wieder benutzt werden können. Sie müssen sich auf einfache 
Art darüber informieren können, wie das Programm aufgebaut 
ist, welcher Lösungsweg gewählt wurde usw. In vielen Fällen 


kann der Programmablaufplan eine spezielle Programmbeschrei- 


bung ersparen. 


2.2 Struktogramme 


Programmablaufpläne sind als Hilfsmittel zur graphischen Dar- 
stellung von Programmen weit verbreitet und auch genormt. Sie 
beschreiben die logische Struktur eines Algorithmus und den 

Ablauf, d.h. die Reihenfolge der Schritte zu programmierender 
Probleme. Die Praxis zeigt jedoch, daß dieses Hilfsmittel den 
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Programmierer dazu verführt, Programme ohne größere Überlegung 
an beliebigen Stellen zu verzweigen und mit Hilfe der Über- 
gangsstellen an beliebigen Orten, oft auf anderen Blättern, 
wieder zusammenzuführen. Dadurch ist ein Programm vielfach 
nicht mehr in einfache, selbständige Blöcke aufteilbar, d.h. 
die Strukturierung des Programms wird erschwert oder sogar un- 
möglich gemacht. Die Strukturierung ist jedoch für umfangrei- 
che Problemstellungen sehr wichtig. Aus diesem Grunde haben 
NASSI und SHNEIDERMANN eine graphische Darstellungsmethode 
entwickelt, die einen "Sprung" von einem Punkt zu einem ande- 
ren Punkt im Programmablaufplan verhindert. Damit wird gleich- 
zeitig sichergestellt, daß nicht nur der "Programmfluß", son- 
dern auch die "Programmstruktur" deutlich wird, d.h. die Be- 
deutung der einzelnen Programmteile für den Gesamtablauf wird 


auf den ersten Blick erkennbar. 


ann Tamm 


Prozeß 
(Aktivitat, Operation) 




















Das Prozeßsinnbild dıent zur Darstellung eınes oder mehrerer 
Befehle wie z.B. 


© Wertzuweisungen (arithmetische Zuordnungsanweiısungen) 

e Ein- und Ausgabeanweisungen (diese besıtzen ın Programm- 
ablaufplanen ein spezielles Sinnbild) 

e Unterprogrammaufrufe (diese besıtzen in Programmablauf- 
planen ebenfalls eın spezielles Sınnbild, das allerdings in 
Bıld 1 nicht aufgenommen wurde). 


Die Form des Prozeßsinnbildes ıst rechteckig. Die Große ist frei 
wahlbar. 



















Das Verzweigungssinnbild dient zur Darstellung bedingter Ver- 
zweigungen mit zwei Alternativen (Ja/Nein-Entscheidung). Das 
Verzweigungssinnbild besteht aus dreı Dreiecken. Das mittlere 
Dreieck (Text) enthalt die Bedingung (Frage), die entweder mıt 
NEIN (F = FALSE) oder JA (T = TRUE) zu beantworten ıst. 
Je nach Beantwortung der Frage wırd das Programm mit einem 
Prozeßsinnbild, das direkt auf das linke Dreieck (F) oder auf das 
rechte Dreieck (T) folgt, fortgesetzt. 


Dıe Große des Verzweigungssinnbildes ist von der jeweiligen An- 
wendung und dessen Erfordernissen abhangig. 


Verzweigung 
(Entscheidung, 
Selektion) 





















Das Wiederholungssinnbild dient zur Darstellung von Schleifen. 
Dıe Wiederholungsbedingung steht links oben im Wiederholungs- 
sinnbild. Hier wird z.B. dıe Zahl der Wiederholungen angegeben 
oder die Bedingung, unter der die Wiederholung zu beenden ist. 


Wiederholung 
(Schleife, Iteration) 





Wiederholungsbed. 


Die zu wiederholenden Anweisungen stehen ım inneren Rechteck 
(Rumpf). Der Rumpf kann aus eıner Struktur beliebiger Ver- 
wicklung bestehen. 


Eine Verschachtelung von Schleifen ist moglıch (weitere Schleı- 
fen im Rumpf!). 


Das Anfang- und Ende Sinnbild dıent zur Darstellung des Be- 
ginns oder des Endes von Programmen. 






Anfang und Ende 





Bi 1d 3: Grund-Symbole für Struktogramme 
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Das Lesen und Zeichnen von Struktogrammen ist schnell erlern- 
bar, da im wesentlichen nur 4 Basis-Symbole verwendet werden, 
die in speziellen Fällen etwas verändert werden. Sie sind zwar 
nicht genormt, jedoch in der Literatur bislang gleichartig 
dargestellt worden. Die 4 Basis-Symbole und zwei häufig benutz- 


te Varianten zeigen Bild 3 und Bild _4. 


Mehrfachverzweigung Das Mehrfachverzweigungssinnbild dient zur Darstellung beding- 

(Schalter) ter Verzweigungen mit mehr als zwei Alternativen. Das obere 
Dreieck des Sinnbildes enthält die Fallabfrage (Bedingung), d.h. 
hier wird angegeben, unter welcher Bedingung zu den einzinen 
Fällen (1, 2,... n) verzweigt wird. 


Schleife mit Dieses Sinnbild dient zur Darstellung von Schleifen, die unter 
Wiederholungsbed. Abbruchbedingung bestimmten Bedingungen abzubrechen sind. Es gleicht dem 


normalen Wiederholungssinnbild mit dem Unterschied, daß im 
Rum Rumpf die Abbruchbedingung aufgenommen ist. 
Abbruchbed. 
EEE 





Bi 1d 4: Häufig benutzte Symbolvarianten für Struktogramme 


2.3 Beispiele 
Beispiel 1 


Es soll die Summe der ganzen geraden Zahlen von 1 bis 20 ge- 
bildet werden. Geben Sie zur Lösung dieses Problems sowohl 


einen Programmablaufplan als auch ein Struktogramm an. 


Zunächst müssen die Anfangswerte gesetzt werden. So wird zu- 
nächst der Inhalt der Speicherzelle, in der die Summen nach 
jedem Schleifendurchlauf abgespeichert werden, Null gesetzt, 
damit Werte, die vorher möglicherweise in der Speicherzelle 


standen, das Ergebnis nicht verfälschen können (Summe = 9). 


Außerdem wird der Anfangswert der Laufvariablen I auf den er- 
sten ganzzahligen Wert, der zur Summe beiträgt, festgelegt 
(I=2). 

Es erfolgt anschließend die erste Rechnung: SUMME = @ + 2. 


Darauf folgt die erste Indexerhöhung, die sich aus I=1I +2 
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Summenwert 
am Anfang 


Wert der Lauf- 
_— variablen am 


Schleifenanfang 





SUMME = 
SUMME + 


_ T Rechenanweisung 





— — 1 Indexerhöhung 


_._ { Schleifenende 


Hinweis: 
Das Symbol & 
bedeutet ‚‚Null‘ 


mit den vorgegebenen Werten zu I = 2 + 2 = 4 ergibt. Durch 
eine Abfrage, ob I < 2® ist, wird der Index I überprüft. Ist 
der Wahrheitswert wahr, wird die zweite Rechnung SUMME = 2 + 4 
ausgeführt usw. Die Schleife wird so lange durchlaufen, bis 
I> 20 wird. Dann wurden alle ganzen geradzahligen Zahlen ein- 


schließlich 20 aufaddiert und die Summe kann gedruckt werden. 


Wiederhole I = 2,4,6,...,2 


SUMME = SUMME + 
Drucke: SUMME 
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Beispiel 2 


Es soll die Ablaufstruktur folgender Lampenschaltung graphisch 


mit Hilfe eines Programmablaufplanes und eines Struktogrammes 


dargestellt werden, die die 


Lampe zum Leuchten bzw. 


B 
nicht zum Leuchten bringt. Zune 
L 
C 
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START 


Schalter A 
geschlossen 




















Schalter B 


geschlossen 
? 





Schalter C 


geschlossen 
? 


Lampe L Lampe L 
leuchtet nicht leuchtet 


BEGIN 


Schalter A geschlossen 

| TT— ? nn 
Schalter B ge- 

schlossen 

? 

Schalter C 

geschl. 

? 


Lampe Lampe 
L L 
leuchtet | leuchtet 
nicht 
















Lampe L leuchtet nicht 











Lampe L 
leuchtet 


END 
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Beispiel 3 


Aus einer bestimmten Anzahl N von positiven Zahlen A(I) soll 
durch einen Suchvorgang die maximale Zahl Amax ermittelt wer- 
den. Die Ablaufstruktur soll mit Hilfe eines Programmablauf- 


planes und eines Struktogrammes dargestellt werden. 


Programmablaufplan: 


Le 


All- Anax> 8 ” 
? 
guu 


Ausgabe: 


© 
® 
8, : 
S 


Amax 


STOP 


Eingabe: N, A (I) 
Amax * ® 


Wiederhole I = 1 bis N 


nn Al) - Amax > 6 eu 


Ausgabe: Amax 
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Ein Spiel - algorithmisch betrachtet 


von Rüdeger Baumann 


Warum Spiele? fragen Sie. 
Ich antworte: um die Kunst 


der Erfindung zu vervollkommnen. 


LEIBNIZ 1716 


Das folgende Beispiel unterbreite ich aus drei Gründen: er- 
stens seines ästhetischen Reizes und der vielfältigen mathema- 
tischen Bezüge wegen; zweitens, weil es die gemeinsame algo- 
rithmische Struktur aller deterministischen Zweipersonenspiele 
zeigt; und drittens, weil es Gelegenheit gibt, Vorzüge der Pro- 


grammiersprache Pascal gegenüber BASIC herauszustellen. 


1 WYTHOFFs NIM 


Ein schon den alten Chinesen bekanntes, von dem Holländer 
W.A. Wythoff 1907 neuentdecktes Nim-ähnliches Spiel lautet 
(in geometrischer Fassung) 
folgendermaßen: 

Gegeben ist ein schach- 
brettartiges Spielfeld 
beliebiger Größe und 


eine Damefigur. 





Spieler A setzt die Dame auf eines der Felder der obersten Zei- 
le oder der am weitesten rechts befindlichen Spalte (in der Ab- 
bildung schraffiert). Die Dame kann wie beim Schachspiel be- 


wegt werden, aber nur in den Richtungen West, Südwest, Süd 
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(siehe die Pfeile). B hat den ersten Zug, dann wird abwechselnd 
gezogen. Wer nicht mehr ziehen kann, weil die Dame im Eck O 


steht, hat verloren. 


Wir stellen sofort fest: Steht die Dame in der Zeile, Spalte 
oder Diagonale, welche die Null O enthält, kann der am Zug 
befindliche Spieler unmittelbar gewinnen, indem er die Dame 
auf das O-Feld zieht. Wir 
markieren die Felder 

dieser Reihen (mit jeweils 
einem Stern *), da wir sie 
beim Setzen der Dame 


vermeiden müssen. 


OD Bd Tr NO So © 


23456789%o0Nm 





Als nächstes nicht markiertes Feld erhalten wir (1:2): es be- 
zeichnet eine Verluststellung für den gerade am Zug befindli- 
chen Spieler. Wie er auch zieht, immer gerät er auf ein mit * 
markiertes Feld und muß 
verlieren. Auf diese Weise 
ermitteln wir, von der End- 
stellung her aufsteigend, 
rekursiv die Verluststel- 
lungen und belegen sie mit 
dem Zeichen O. 


do rer no SS 9 © 


012345678990 
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Es handelt sich um eine Art Siebverfahren (analog zum Sieb des 


Eratosthenes), das sich folgendermaßen beschreiben läßt: 


WYTHOFF-SIEB 


FÜR x VON ® BIS n WIEDERHOLE 
FÜR y VON ® BIS n WIEDERHOLE 
WENN brett(x,y) = ® DANN 
gib position x,y aus; (* Verluststellung *) 


markiere alle von x,y ausgehenden waagerechten, 
senkrechten und diagonalen reihen 
ENDE-WENN 
ENDE-WIEDERHOLE 
ENDE-WIEDERHOLE 





Dies ist der Kern einer vom Computer anzuwendenden Strategie: 
er braucht die Dame lediglich auf die nächstliegende so ermit- 


telte Verluststellung zu ziehen (sofern sie nicht schon auf 
einer steht; in diesem Fall macht der Computer einen Verlegen- 


heitszug und hofft auf einen Fehler des Gegners). 


2 SPIELPROGRAMM 


Jedes Wettkampfspiel zwischen zwei Personen (Mensch und Compu- 
ter) läßt sich gliedern, wie es das folgende Baumdiagramm be- 


schreibt (siehe Schrage [14]): 


Spiel 
Spiel- Start- Partie Gewinn- 
beschreibung angaben (Zugfolge) entscheid 
Anzeige zug des zug des Ende- 
des Spielstands Menschen Computers prüfung 
Eingabe Prüfung Bestimmung Zug- Bestimmung Ausgabe 
des Zugs des Zugs des neuen ermitt- des neuen des Zugs 
auf Zuläs- Spiel- lung Spiel- 
sigkeit stands stands 
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Die Besonderheiten unseres Spiels kommen erst bei der Ermitt- 
lung des Computerzugs zur Geltung. Bevor wir darauf eingehen, 
wird - gemäß dem Baumdiagramm - die Grobstruktur des Spielpro- 


gramms entwickelt; es lautet in der Programmiersprache Pascal: 


PROGRAM treib-die-dame; 


(* nim-ähnliches strategiespiel, das von den alten chine- 
sen erfunden, und von w.a. wythoff neuentdeckt wurde *) 


CONST 


seitenlänge 50; 


TYPE 


bereich 
feldtyp 
spielbrettyp 
stellungstyp 
spielertyp 


Ö..seitenlänge; 

(leer, markiert, dame); 
ARRAY[bereich, bereich] OF feldtyp; 
RECORD x,y : bereich END; 

(mensch, computer); 


VAR 


brett spielbrettyp; 

position stellungstyp; 

amzug spielertyp; (* kennzeichnet den am Zug 
befindlichen Spieler *) 


PROCEDURE beschreibung; 
(* gibt die spielanleitung und benutzerhinweise *) 
BEGIN (* wird hier nicht ausgeführt *) END; 


PROCEDURE startangaben; 
(* legt anziehenden spieler und startposition fest *) 
BEGIN (* wird hier nicht ausgeführt *) END; 


PROCEDURE partie; 
(* das eigentliche spiel, die zugfolge *) 
BEGIN (* wird weiter unten ausgeführt *) END; 


PROCEDURE gewinnentscheid; 
(* ermittelt den gewinner der partie *) 
BEGIN (* wird weiter unten ausgeführt *) END; 


BEGIN (* Hauptprogramm *) 
beschreibung; 
startangaben; 
partie; 


gewinnentscheid 


END. 
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In der Prozedur "Startangaben' wird der anziehende Spieler be- 
stimmt, und es wird die Anfangsposition der Dame durch jenen 
festgelegt. Wegen der Platzbeschränkung in diesem Band muß ich 
auf den Abdruck der Prozedur verzichten; interessierten Lesern 
kann ich einen Ausdruck des vollständigen Pascal- (oder auch 
BASIC-)Programms zusenden. Die Prozedur '"Gewinnentscheid' ist 


ganz einfach: 


PROCEDURE gewinnentscheid; 
(* ermittelt den gewinner der partie *) 
BEGIN 


IF amzug = mensch THEN 
writeln('Sie haben verloren.') 


ELSE 
writeln('Sie haben gewonnen!') 


END; (* of gewinnentscheid *) 





Interessant ist allein die Prozedur 'Partie'; ihre Grobstruk- 


tur lautet wie folgt: 


PROCEDURE partie; 
(* das eigentliche spiel, die zugfolge *) 


VAR 
links, unten : bereich; 


PROCEDURE spielstandsausgabe; 
(* gibt den aktuellen spielstand aus *) | 
BEGIN ( * wird im folgenden nicht ausgeführt *) END; 


PROCEDURE neuer-spielstand; 
(* der neue spielstand wird berechnet *) 
BEGIN (* im folgenden nicht ausgeführt, da sehr ein- 
fach *) END; 


PROCEDURE der-mensch-zieht; 
(* der spielzug des menschen wird eingelesen und aus- 
geführt *) 
BEGIN (* wird im folgenden dargestellt *) END; 


PROCEDURE der-computer-zieht; 
(* der spielzug des computers wird ermittelt und aus- 
geführt *) 
BEGIN (* wird im folgenden genauer dargestellt *) END; 


FUNCTION partie-ende : boolean; 
BEGIN 
partie-ende := (position.x = ®) AND (position.y = $) 
END; (* of partieende *) 
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BEGIN (* partie *) 
REPEAT 


spielstandsausgabe; 


IF amzug = mensch THEN BEGIN 
der-mensch-zieht; 


amzug := computer END 
ELSE BEGIN 
der-computer-zieht; 
amzug := mensch 
END (of else*) 


UNTIL partie-ende 
END; (*of partie*) 





Die Prozedur "der Mensch zieht" lautet wie folgt: 


PROCEDURE der-mensch-zieht; 
(* der spielzug des Menschen wird eingelesen und aus- 
geführt *) 


FUNCTION in-ordnung : boaolean; 
(* prüft, ob der spielzug legal ist *) 
BEGIN (* wird nicht näher ausgeführt *) END; 


PROCEDURE zugeingabe; 


(* der zug des menschen wird eingelsen *) 
BEGIN (* nicht ausgeführt, da einfach *) END; 


BEGIN 


REPEAT zugeingabe UNTIL in-ordnung; 
neuer-spielstand; 


END; (* of der mensch zieht *) 
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Analog dazu die Prozedur "der Computer zieht": 









PROCEDURE der-computer zieht; 
(* der spielzug des computers wird ermittelt und aus- 
geführt *) 









VAR 
x,y : bereich; 










PROCEDURE zugausgabe; 
(* der ermittelte spielzug wird ausgegeben *) 
BEGIN (* hier nicht ausgeführt, da klar *) END; 










PROCEDURE zugermittlung; 
(* der computer .ermittelt seinen besten zug *) 
BEGIN (* wird unten näher beschrieben *) END; 








BEGIN 


zugermittlung; 
neuer-spielstand; 
zugausgabe 








END; (* of der computer zieht *) 






Die Spielstrategie steckt in der Prozedur 'Zugermittlung': 


PROCEDURE zugermittlung; 
(* der computer ermittelt seinen besten zug *) 


VAR 
i : integer; 


PROCEDURE waagerechtmarkieren" 
(* eine zeile wird markiert und ggf. die dame entdeckt *) 
BEGIN (* wird unten näher erläutert *) END; 


PROCEDURE senkrechtmarkieren; 
(* eine zeile wird märkiert und ggf. die dame entdeckt *) 
BEGIN (* analog zu waagerechtmarkieren *) END; 


PROCEDURE diagonalmarkieren; 
(* eine diagonale wird markiert und ggf. die Dame ent- 
deckt *) 
BEGIN (* analog zu senkrechtmarkieren *) END; 





BEGIN 
FOR x := Ö TO seitenlänge DO 
FOR y := ® TO seitenlänge DO 
brett[x,y] := leer; 
brett[x,y] := dame; 
links := ®; unten := $; 
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FOR x := ® TO seitenlänge DO 
FOR y := Ö TO seitenlänge DO BEGIN 
IF brett[x,y] = leer THEN BEGIN 
writeln(x,' ',y); (* verluststellung *) 
waagerechtmarkieren; 
senkrechtmarkieren; 
diagonalmarkieren 
END; 
IF (links = ®) AND (unten = ®) THEN 
links := 1 (* Verlegenheitszug *) 
END (* of for *) 
END; (* of zugermittlung *) 





Es bleibt die Prozedur "waagerecht markieren", die die von 
einer Verluststellung ausgehende Zeile markiert; findet sie 
auf einem Feld die Dame, so zieht sie diese waagerecht auf die 
Verluststellung, von der sie ausging und springt aus der Pro- 


zedur 'Zugermittlung' heraus. 


PROCEDURE waagerechtmarkieren; 


(* eine zeile wird markiert und ggf. die dame entdeckt *) 


WHILE x+i <= seitenlänge DO BEGIN 
IF brett[x+i,y] = dame THEN BEGIN 


links := i; unten := 6; 


exit (zugermittlung) 
END; (* of then *) 
brettlx+ti,y] := markiert; 
ji := i+1 
END (* of while *) 


END; (* of waagerechtmarkieren *) 





3 VERGLEICH BASIC-Pascal 


Beim folgenden BASIC-Programm wurde versucht, genau die eben 
beschriebene Struktur (Gliederung in Hauptprogramm und Proze- 
duren) nachzuahmen. Trotzdem wird längst nicht der Grad an Les- 
barkeit und Verständlichkeit wie beim Pascal-Programm erreicht; 


u.a. aus folgenden Gründen: 


- Prozeduren lassen sich in BASIC nicht benennen und unter ih- 


rem Namen aufrufen, 
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- Variablen haben in BASIC nur zwei signifikante Zeichen, 
- die Definition von Datentypen, wie z.B. unser 
"feldtyp = (leer, besetzt, dame)" ist in BASIC nicht möglich. 


Ein großer Vorteil von BASIC gegenüber Pascal soll indessen 
nicht verschwiegen, sondern muß gebührend herausgestellt wer- 
den: BASIC ermöglicht schon zwölfjährigen Schülern, komplexe 
Spielprogramme (wie z.B. das vorliegende) zu entwickeln; in 
Pascal würden sie dagegen bereits an den Deklarationen schei- 
tern. Ihr Arbeitsstil ist intuitiv und experimentierend, durch 
Versuch und Irrtum tasten sie sich an das ihnen vorschwebende 
Ziel heran: dergleichen ist in Pascal nicht möglich. 


PRINT CHR$ (147) 
PRINT " TREIB DIE DAME IN DIE ECKE 


PRINT!" | -4a44444 


1 Ö OÖ (®) 
1010 
1011 
1012 : 
1020 REM NIM-AEHNLICHES STRATEGIESPIEL, DAS VON DEN ALTEN CHINESEN 
1021 REM ERFUNDEN UND VON W.A.WYTHOFF (1907) NEU ENTDECKT WURDE 
1022 3 

1099 : 
1100 REM KK HAUPFTPROGRAMM RRRRRORRRRRRRRORIORIDDRIRRIDRRRDRRIRRRRKIOKEK 
1101 : 

1200 


GOSUR 2000 REM BESCHREIBUNG 


13200 2 GOSUR 3000 : REM STARTANGABEN 

1400 : GÜSUR 4000 = REM PARTIE 

1500 » GOSUR 5000 s REM GEWINNENTSCHEID 

1900 : END 

1901 ; 

1990 REM ENDE DES HAUPTPROGRAMMS ORRRRRIRRRIRRRIIRRRRRRIORRDIGOIIKEK 
1998 : 

1999 ı 

2000 REM ### PROZEDUR °’BESCHREIRBUNG’ ####HHHHH HH 
2001 : 

2100 2 REM HIER SOLLTEN DIE SPIELREGELN UND BENUTZERHINWEISE STEHEN 
2101 : 

2900 » RETURN 

2901 : 

2990 REM ENDE DER PROZEDUR "BESCHREIBUNG” #HHHH#HHH tt 
2999 : 


S000 REM ### FRÜZEDUR "STARTANGABEN’ HEHE 
sool 


050 = LETN = 20 : REM SEITENLAENGE DES BRETTS 

=100 : DIM B(N,N) : REM SPIELERETT 

190 : 

200 2 FRINT "GEBEN SIE DIE ANFANGSFOSITION EIN! 

3201 3: 

205 s PRINT 

3210 » INPUT "X-KDORDINATE "3 X 

S220 3 INPUT "Y-EKDORDINATE "3 Y 

Sal: 

230 2 IFX => N ORY > N THEN PRINT "ZU GROSS!"s GOTO 7210 
S240 : LET B(X,Y) =5 

00 3 LET SF$ = "COMPUTER" 

S>501 : 

3900 s RETURN 

S>901 3: 

3990 REM ENDE DER PROZEDUR °STARTANGABEN' HH He 
3999 : 
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4000 
4001 

4100 
4101 

4110 
4120 
4121 

4130 
4131 

4135 
4140 
4145 
4149 
4150 
4151 

4158 
4155 
4160 
4165 
4170 
4171 

4175 
4180 
4185 
4190 
4191 
4195 
4199 
4200 
4201 
4210 
4220 
4221 
4290 
4291 
4295 
4299 
4300 
4301 
4310 
4320 
4330 
4340 
4350 
4360 
4361 
4370 
4380 
4385 
4389 
4390 
4391 
4395 
4399 


4400 
4401 
4405 
4406 
44067 
4409 
4410 
4420 
4430 
4435 
4440 
4450 
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REM ### PROZEDUR "’PARTIE’ SEHHH HEHE 


REM === PROZEDURRUMFF ’FARTIE” =====2=: =2=2=-=m 22 2- = 2=2mese=msmeuememecmenm 


GOSUB 4200 =: REM SPIELSTANDSAUSGABE 
IF SP$ = "COMFUTER" THEN 4150 


REM --- DER MENSCH IST AM ZUG 
GOSUB 4500 =: REM ZUG DES MENSCHEN 
LET SP$ = "COMPUTER" 
GOTDO 4170 : REM PARTIE ZU ENDE? 
REM --- DER COMPUTER IST AM ZUG 
PRINT s PRINT "LASSEN SIE MICH UEBERLEGEN ... 
GOSUB 4400 : REM ZUG DES COMPUTERS 
LET SP$ = "MENSCH" 
REM --- PRUEFUNG AUF PARTIEENDE 


GOSUB 4900 
IF E$ = "WEITER" THEN 4110 


RETURN 


REM ENDE DES PROZEDURRUMFPFS "PARTIE’ =========2=2=2==2=2=2==22220002= 


REM +++ UNTERFROZEDUR ”SPIELSTANDSAUSGABE? +++++++++++tttt+t++++++ 


PRINT 
PRINT "SPIELSTAND: "; X; Y 


RETURN 


REM ENDE DER UNTERPROZEDUR °’SPIELSTANDSAUSGABE’ ++++++++++++++++ 


REM +++ UNTERPROZEDUR ’ZUG DES MENSCHEN’ ++++4++ +++ tt tt tt tt HH 


PRINT 

FRINT "SIE SIND DRAN! 

INPUT "WIEVIEL NACH LINKS "5 L 

INPUT "WIEVIEL NACH UNTEN "3 U 

IF LXU <> OO AND L <> U THEN PRINT "FALSCHER ZUG!"s GOTO 4330 
IFL = 0 AND U = © THEN PRINT "SIE MUESSEN ZIEHEN!": GOTO 4330 


LET X = X-L : LET Y = Y-U 
IFX< 0 0O0RY< 0 THEN 4330 
LET B(X,Y) =5 


RETURN 


REM ENDE DER UNTERPROZEDUR ”’ZUG DES MENSCHEN’ +++++++++++++++++ 


REM +++ UNTERPRÜOZEDUR ?’COMPUTERZUG’? ++++++++ tt tt tt Hr tr tr tr tr tr 


FORR=0TON: FORS=0TON : LET B(R,S) = 0 2 NEXT S;R 
LET B(X,Y) =5 
LETL=0: LETU=0© 


FORR = 0 TON 
FOR S=0TON 

IF B(R,S) = 

PRINT R,S : 


1 THEN 4400 
REM AUSGARE DER VERLUSTSTELLUNG 


LET I=0 
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4460 : LET I = I+i1 : IF R+I > N THEN 4500 

4470 : IF B(R+I1,5) <> 5 THEN LET B(R+1,5) = 1 : GOTO 4460 

4490 LET L = I: LET U = 0 : GoOTO 4700 

4490 : 

4500 5; LET I=0 

4510 ; LET I = I+ti » IF S+I > N THEN 4550 

4520 3; IF B{R,S+I) <> 5 THEN LET B(R,S+I) = 1 : GOTO 4510 

45230 2 LETL = 0: LET U= 1 : GOTO 4700 

4540 :; 

4530 5 LET I =0 

4560 : LET I = I+1 : IF R+I > N OR S+I > N THEN 4400 

4570 : IF B(R+I,S+I) <> 5 THEN LET B(R+I1,S+I) = 1 : GOTO 4560 
4580 3 LETL = I: LET U= I s G0OTO 4700 

4590 3 

4600 3 NEXT S 

4610 : NEXT R 

4699 : 

4700 : IFL = 0O AND U = © THEN LET L = 1 : REM VERLEGENHEITSZUG 

4710 3 PRINT "ICH ZIEHE"; L;3 "NACH LINKS UND"; Us "NACH UNTEN. 

4720 3: LET X = X-L : LET Y = Y-U 

4720 3 

4790 : RETURN 

4791 3: 

4795 : REM ENDE DER UNTERFRÜZEDUR "COMPUTERZUG’? +r++++++ttrr tr tr tr tr Hr rt + 
4799 3 

4900 s REM +++ UNTERFROZEDUR ’PARTIEENDE” +++++++++++t++ttttttttrtr tt rt HH 
4901 : 

4910 ı LET E$ = "WEITER" 

4920 :; IFX = 0 AND Y = 0 THEN LET E$ = "ENDE" 

4930 :; 

4990 3 RETURN 

4991 : 

4995 : REM ENDE DER UNTERPRÜZEDUR "PARTIEENDE’ +++++++++++++++++t++++++ 
4999 : 

SO00 REM ##4# FPROZEDUR ”’GEWINNENTSCHEID’ #EHHHH HH HH HH HH HH HH HEHE HN 
soo0l : 

3100 s IF SP$ = "MENSCH" THEN PRINT "SIE HABEN VERLOREN. 

5200 : IF SF$ = "COMPUTER" THEN FRINT "SIE HABEN GEWONNEN! 

2:00 3: 

5900 3 RETURN 

s9a0l ; 

990 REM ENDE DER PROZEDUR ’GEWINNENTSCHEID’ #HHHHH HH 
3999 3 


4 MATHEMATISCHE ENTDECKUNGEN 


Beim Spiel gegen den Computer produziert dieser auf seiner Su- 
che nach einem optimalen Zug die Folge der Zahlenpaare x(n), 


y(n), welche die Verluststellungen kennzeichnen: 










1 213] 4| 5| 6 8| 9j 10] 11 | 12| 13| 14 16 
JEIEICIEEIN ZIE 
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An ihr lassen sich interessante mathematische Entdeckungen ma- 


chen. Wir stellen zunächst fest: 
(1) x(n) + y(n)=n 


(2) x(n) ist die kleinste der bisher noch nicht aufgetretenen 
zahlen. 


Hieraus läßt sich die Folge der Verluststellungen ohne die auf- 
wendige Suche in der Prozedur "Zugermittlung' rein numerisch 
(rekursiv) gewinnen (und zu Beginn eines Spiels speichern). 


Dies ist wieder ein Beispiel für den Sachverhalt: 


- Mit etwas Mathematik kann man sich aufwendiges Programmieren 
ersparen. 


Sein Gegenstück lautet: 


- Mit etwas Programmieren kann man (eventuell schwierige) Ma- 
thematik umgehen; 


beide zusammen kennzeichnen das ambivalente Verhältnis von Com- 
puter und Mathematik, welches zunehmend den Unterricht zu be- 
herrschen beginnt. - 

Dem Kenner fällt sofort auf, daß in obiger Tabelle die Fibo- 
nacci-Folge versteckt ist; damit ist der Weg zum goldenen 
Schnitt nicht weit (vgl. [4]). Die Folgen x(n), y(n) sind in 


der Form 
(3) x(n) = int (n-y) 
(4) y(n) = int (n-o%) mit 9= 1+2 


explizit darstellbar (die Programmieraufgabe ist jetzt schon 
fast trivial!). Hieraus läßt sich eine direkte Kennzeichnung 
der Verluststellungen entwickeln, die an die berühmte Dualdar- 
stellung der Verluststellungen im Nim-Spiel erinnert. 

Aus Platzgründen kann ich diese interessanten - etwa für eine 
Mathematik- oder Informatikarbeitsgemeinschaft bzw. für den 
sogenannten problemorientierten Unterricht (vgl. [8]) geeigne- 
ten - Möglichkeiten für mathematische Entdeckungen nicht näher 
eingehen; genaueres findet der Leser in [6], und weitere Lite- 
ratur in [3, S. 80]. 
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Zwei Knobeleien 


von Rüdeger Baumann 


Die folgenden Beispiele sollen zeigen, 


- daß Pascal wohlstrukturiertes Programmieren unterstützt, 
während BASIC zum sogenannten Spaghetti-Code verführt, 
- daß Pascal zur Modellierung realer Situationen geeignete Da- 


tenstrukturen besitzt, während BASIC hieran arm ist. 


1. n-DAMEN-PROBLEME 


Ein klassisches Problem der Unterhaltungsmathematik lautet; 


Auf einem Schachbrett sollen acht Damen so plaziert wer- 





den, daß sie sich gegenseitig nicht schlagen können. 





Selbst der große C.F. Gauß hat sich mit dieser Aufgabe beschäf- 
tigt; und natürlich interessierte er sich nicht nur für eine, 

sondern für alle Lösungen, mehr noch: er suchte ein systemati- 
sches Verfahren zur Erzeugung aller Lösungen. Es läßt sich wie 


folgt beschreiben: 


"Man setzt von links nach rechts fortschreitend in jede Spalte 
des Schachbretts eine Dame, und zwar jedesmal an die tiefst- 
mögliche Stelle. Wenn nun der Augenblick kommt, wo man keine 
Dame in einer Spalte aufstellen kann, so erhöht man den Platz 
der Dame in der vorhergehenden Spalte so lange, bis man fort- 


fahren kann." 


Die Stellung der Damen auf dem Schachbrett kennzeichnen wir 
durch Zahlen D(i); dabei ist i die Nummer der jeweiligen Spal- 
te, in der die Dame steht (i = 1,2,...,n). Ein Flußdiagramm 


zum Gaußschen Verfahren ist in B i 1 d 1 wiedergegeben. 
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Lösung drucken 


Di):=Di)+1 


Bild 1 1: Flußdiagramm zum n-Damen-Problem 
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Das zugehörige BASIC-Programm lautet: 


100 INPUT „Wieviele Damen’; n 








110 DIM Din) 
120 LETi=O 
| 130 LETi=si+i 
| 140 LET Di) = 1: REM erste Dame aufs Brett 
150 IFi=1THEN 190 
160 FORk= 1 TOi 1: REM Bedrohung prufen 
| 170 IF Di) = Dik) OR ABS(DÜi) Dik))=i KkTHEN 230 
| 180 NEXTk 
190 IFi<n THEN 130 
200 FORk=1TOn PRINT D(k);: NEXTk 
210 PRINT : REM Lösung gefunden und gedruckt 
220 LETi=i 1: REM Dame vom Brett nehmen 
| 230 LET D(i) = Dii) + 1: REM Dame ruckt vor 
240 IF D(Üi) < =n THEN 150 
250 LET i = ı - 1: REM Dame vom Brett 
260 IFi > 0 THEN 230 
zwei Damen in Spalte i und k bedrohen sich, wenn D(i) = D(k) 


ist (dann stehen sie in der gleichen Zeile) oder wenn 
|ID(i)-D(k)| = Ji-k| (gleiche Diagonale): dies wird in Zeile 
170 überprüft. 


Das Programm ist erstaunlich kurz; dies verdankt es vor allem 
der in BASIC gegebenen Möglichkeit, an jede beliebige Stelle 
innerhalb eines Programms zu springen. Doch birgt sie Gefahren, 
denn das Springen (Anweisung GOTO), vor allem das Springen in 
Schleifen (wie es das Flußdiagramm anschaulich zeigt), macht 
Programme fehleranfällig und erschwert die Fehlersuche. 
Schlimmer noch: Programme dieser Art werden ab einem gewissen 
Umfang völlig unübersichtlich und ein Beweis der Korrektheit 
(hier: Beweis der Tatsache, daß mit dem Programm alle Lösungen 


des Acht-Damen-Problems gefunden werden) ist fast unmöglich. 


Wollen wir das Gaußsche Verfahren in der Programmiersprache 
Pascal darstellen, müssen wir die 'Spaghetti-Fäden' zunächst 


entwirren, d.h. dem Algorithmus eine saubere Struktur geben. 


Sie läßt sich am besten in einem Struktogramm ausdrücken (zu 
dem obigen Flußdiagramm dagegen gibt es kein Struktogramm!). 
Dies ist in Bi ld 2 angegeben. 


Als Preis für die Wohlstrukturiertheit müssen wir die doppel- 
te Abfrage "D(i) > Spaltenzahl" bzw. "D(i) < Spaltenzahl" zah- 


len; sie kam im BASIC-Programm nur einmal vor (Zeile 240). Die 
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i:=1  (* Spaltenindex wird initialisiert*) 
D(i) :=0 (* Dame Nr. 1 steht noch außerhalb des Bretts*) 


i=0 (* alle Lösungen sind gefunden*) 







Di) :=Dii) +1 (* Dame rückt einen Platz vor*) 


Dame nicht bedroht 


7 D(i) > Spaltenzahl (* Dame außerhalb des Bretts*) 


Di) <S Spaltenzahl . 
R i < Spaltenzahl 


jet 
(* neue Dame Lösung drucken 
aufs Brett*) 


om:=0 _|tdam 


(* Dame wird vom Brett genommen*) 














Bi 1d 2: Struktogramm zum n-Damen-Problem 
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Blöcke "Dame nicht bedroht" und "Lösung drucken" legen wir in 


Unterprogramme (Prozeduren) : dies fördert die Übersichtlich- 


keit. 
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Das Struktogramm läßt sich unmittelbar in Pascal übersetzen: 


PROGRAM achtdamen; 


CONST 
spaltenzahl = 8; 
VAR 
brett : ARRAY |1.. spaltenzahl | OF integer; 
index: integer; (* Spaltenindex der Damen *) 
FUNCTION bedroht (1: integer) : boolean; 
VAR 
k : integer; 
BEGIN 
bedroht := false; 
FORk :=1TOi-1DO 
IF (brett ji] = brett |k ]) OR (abs (brett [i] brett Ik) = i 
THEN BEGIN 
bedroht := true; exıt (bedroht) 
END (* ofthen *) 
END; (* of bedroht *) 
PROCEDURE ausgabe; 
VAR 
k : ınteger; 
BEGIN 
FOR k := 1 TO spaltenzahl DO write (brett |k] : 4); 
writeln 
END; (* of ausgabe *) 
BEGIN (* Hauptprogramm *) 
ındex .= 1; brett |index] = O; 
WHILE index > 0 DO BEGIN 
REPEAT 
brett |index] := brett [index] + 1 . 
UNTIL NOT bedroht (index) OR (brett |ındex] > spaltenzahl); 
IF brett [index] < = spaltenzahl 
THEN 
IF ındex < spaltenzahl 
THEN BEGIN 
index = ındex + 1; (* neue Dame *) 
brett [index | = O END (* of then *) 
ELSE BEGIN 
ausgabe; (* Lösung gefunden *) 
ındex := ındex 1 (* Dame vom Brett *) 
END (* ofelse *) 
ELSE ındex ‘= ındex 1 (” Dame vom Brett *) 
END (* of while *) 
END 





Mit diesen Programm lassen sich die 92 Lösungen des Acht-Da- 
men-Problems, angefangen bei 15863724, bis 84136275 erzeugen; 
dabei läuft das Pascal-Programm auf einer Microengine etwa 


100mal so schnell wie das BASIC-Programm auf einem CBM. 


Beachtet man, daß sich aus einer Lösung durch Drehen und Spie- 
geln weitere erzeugen lassen, reicht es, den Computer nur bis 
D(1) = 3 rechnen zu lassen: dies ermöglicht eine drastische 
Verkürzung der Rechenzeit. Übrigens kann das Pascal-Programm 
auch rekursiv formuliert werden; eine Möglichkeit, die es in 


BASIC nicht gibt. 
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2 JOSEFS-SPIEL 


Eine besondere Stärke von Pascal ist das reiche Angebot an 
höheren Datenstrukturen: es gibt - wie in BASIC - das Feld 
(array) ; darüber hinaus aber Menge (set), Verbund (record), 
sequentielle Datei (file) sowie den Zeigertyp (pointer). Letz- 
terer ist besonders interessant, weil er die Möglichkeit bie- 
tet, Variablen dynamisch zu kreieren. Als Exempel hierfür die- 
nen das sogenannte Josefs-Spiel [1], ebenfalls ein altes Pro- 
blem der Unterhaltungsmathematik, welches in allgemeiner Fas- 


sung wie folgt lautet: 


n Personen stehen im Kreis, jeder k-te wird ausgeschie- 






den, wobei sich der Kreis sofort wieder schließt. Ge- 


sucht ist die Reihenfolge der Ausgeschiedenen. 


In BASIC lösen wir es so: Wir denken uns die Kriegerschar (in 
der alten Geschichte, die dem Problem zugrundeliegt, handelt 

es sich um Soldaten) als verkettete Liste, indem wir in ein Feld 
mit n Plätzen die Nummer des jeweils rechts stehenden Krie- 


gers, also des Nachfolgers in der Zählrichtung, schreiben: 


k(1) k(2) K(3) ... k(n-1) k(n) 


Das heißt z.B.: Krieger 3 ist Nachfolger von Krieger 2, und 
Krieger 1 ist Nachfolger von Krieger n: der Ring ist geschlos- 


sen. 


Das Ausscheiden eines Kriegers wird nun durch die Anweisung 
k(z) := k(k(z)) realisiert: auf Platz z des Feldes steht jetzt 
nicht mehr der alte Nachfolger k(z), sondern dessen Nachfolger 
k(k(z)), d.h. kK(z) wurde ausgeschieden. Er wird bei jedem 
künftigen Durchlaufen des Feldes übersprungen. Das Ausscheiden 
wird so lange fortgeführt, bis nur noch ein Krieger übrig ist: 
er ist dann sein eigener Nachfolger. Die Abbruchbedingung lau- 
tet also k(z) = z. Damit läßt sich das BASIC-Programm unmittel- 


bar hinschreiben: 
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PRINT“ Josefs-Spiel’’ 
PRINT“ “ 





.unasuunnnenennne 


110 REM Demonstriert verkettete Listen 
111 : 
120 REM * Eingabe * 









121 : 

130 INPUT „Wieviel Krieger‘; n 

140 INPUT ‚‚Der wievielte soll jeweils ausgezählt werden’; s 
190 : 

200 REM * Aufstellen im Kreis * 

201 : 





210 DIM k(n) 

220 FOR z=1TOn 1:LETk(z)=z+1:NEXTz 

230 LET kn) = 1 : REM Schließen des Rings 

290 : 

300 REM * Auszählen * 

301 : 

310 LET z=n: REM Ruückstellen des Zeigers z 

320 FORi=1T0Os -1:LETz=k(z):NEXTi 

330 PRINT k(z); 

340 LET k(z) = k(k(z)) : REM Ausscheiden der Nr. z 

350 : 

360 IFkKiz) <> 2 THEN 320 : REM Noch nicht alle ausgezählt 
370 PRINT z : REM Dies ist der letzte 

END 













Dialog: 










Wieviel Krieger? 40 
Der wıevielte soll jeweils ausgezahlt werden? 7 
ı 14 21 28 35 2 ı0 18 26 34 3 12 
22 33 40 1 23 3 5 ı7 30 4 ı9 36 
9 27 625 8 32 16 ı 38 37 39 15 
29 13 20 24 





In Pascal könnte man dieses Programm natürlich mittels des Da- 
tentyps ARRAY nachvollziehen; eine wesentlich elegantere Lö- 
sung bieten aber Variablen vom Zeigertyp. (Aus Platzgründen 
ist es hier leider nicht möglich, den Begriff des Zeigertyps 


ausführlich zu entwickeln; man vergleiche dazu [2] oder [1].) 


Das zugehörige Pascal-Programm ist auf der nächsten Seite 


abgedruckt. 


Dieses Programm hat vor dem BASIC-Programm den Vorteil, daß 
die Anzahl der aufzustellenden Personen vorher nicht bekannt 
zu sein braucht. Weiterhin können wir die Namen der Personen 
(nicht nur ihre Nummern) eingeben: dies wird durch die Benut- 
zung des Datentyps Verbund (RECORD) ermöglicht. Er ist beson- 
ders für die Dateibearbeitung wichtig, existiert aber in BASIC 


leider nicht. 


32 


Zwei Knobeleien 





Pascal-Programm "Josephsspiel": 


———— 


PROGRAM josefsspiel; 





TYPE 
zeıger = } soldat; 
soldat = RECORD 
name : string; 
nachfolger : zeıger 
END; (* of soldat *) 


VAR 
erster, krieger : zeıger; 
schrittweite : ınteger; 


PROCEDURE aufstellenımkreis; 


PROCEDURE erzeuge (VAR k : zeiger); 
BEGIN 
new (k}; 
writeln (’Name?’); 
read (k t name) 
END; (* of erzeuge *) 
BEGIN (* of aufstellenımkreis *) 
erzeuge (erster); 
krıeger := erster; 
REPEAT 
erzeuge (krıeger ! . nachfolger); 
krieger := krıeger ! . nachfolger 
UNTIL eof; 
krieger t . nachfolger := erster 
END; (* of aufstellenimkreis *) 
PROCEDURE auszahlen (schritt : ınteger); 
VAR 
I: Integer; 
BEGIN 
REPEAT 
FOR ı := 1 TO schritt 1DO 
krıeger := krıeger t . nachfolger; 


write (krıeger ? . nachfolger ? . name : 6); 


krieger t . nachfolger := krieger t . nachfolger ! . nachfolger 


UNTIL krieger t . nachfolger = krieger; 
write {krieger t . name 6) 
END; (* of auszahlen *) 
BEGIN (* Hauptprogramm *) 
writeln; 
wrıteln(’ Josefs-Spiel’); 
writein(’ ; 
writeln; wrıteln; 
wrıteln(’Geben Sie die Namen eın'’); 
aufstellenimkreis; 


wrıteln(’Der wıevielte soll jeweils ausgezählt werden?‘); 


read(schrittweite); 
auszahlen(schrittweite) 
END. 
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10 Beispiele 


von Dietmar Herrmann 


1  QUERSUMME 


Ganzzahlige Division durch 1@ und Restbildung trennt die 


Einerziffer einer Dezimalzahl ab. 


2.B. gilt: 
7249 mod 1® = 9 
7249 div 18 = 724. 


Setzt man das Verfahren mit dem Quotienten fort, so kann man 
schrittweise alle Ziffern abtrennen und die Quersumme berech- 
nen. 

Dieses Vorgehen kann direkt in Pascal programmiert werden: 


quersumme := ß; 
repeat 
quersumme := quersumme + n mod 19; 
n := n div 19 
untiln = @; 
Für die Funktion quersumme muß jedoch noch eine Variable 
quers eingeführt werden, da sich sonst die Funktion selbst 
aufruft und somit rekursiv berechnet wird. 
Zur Ausführung der MOD- und DIV-Funktionen muß die Zahl n als 
ganzzahlig vereinbart werden. Dies schränkt die Größe der Zahl 


ein; bei dem hier benutzten Compiler gilt nämlich 
MAXINT = 32767. 


In BASIC können die Funktionen MOD und DIV mit Hilfe der INT- 


Funktion dargestellt werden: 


a mod b entspricht a - int (a/b)*b 
a div b entspricht int(a/b). 


Da man aber in BASIC mit Hilfe der MID$-Funktion alphanumeri- 
sche Variablen in einzelne Zeichen verwandeln kann, liegt ein 


anderes Vorgehen nahe: 
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Quersumme 





Die Zahl Z wird durch 
2$ = STR$(Z) 


in eine Ziffernkette verwandelt. Die Quersumme kann dann wie 


folgt berechnet werden: 


Qa=% 

FOR I=1 TO LEN (Z$) 

Q =0Q + VAL(MID$(2Z$,I,1)) 
NEXT I 


Die VAL-Funktion wandelt die Ziffern wieder in Zahlen um. Da 
die meisten BASIC-Interpreter alphanumerische Variablen bis 
zu 255 Zeichen Zulassen, kann man damit auch die Quersumme 


sehr großer Zahlen berechnen. 
NIE 


100 perceram er summelinput.cautput); 

‚118 var zahl : inteser; 

120 ı# # 
120 functicar Auer zunmmeinsinteser > sintesger ; 
142 var erzs:sinteger; 

150 kegirı 


1e0 wierz:s=ü; 

170 repest 

1230 quer z;=werstn mad 10; 

134 ns=n div 1 

zag until rn=d; 

21a yuer Zufme sZzawuierz 

2zu endzi* werzunmme #) 

230 u#% #5 

240 Kkesir 

25a uriteins "Berechnung der Quersumme” >; 
250 write ln; , 

eva uritelni"gik natver liche Zahl ein!” 
zeug readlzahl>;: 

228 uritelne "die Quersumme von” „zahl:sd. "=" „Wwerzsummelzahl::2) 


zaa erd. 
In 
1a REM BUERSUMME EINER ZAHL 

11a PRIHTCHREN 14T) 

1zu : 

120 IHMPUT"EIHNGABE DER ZAHL" :Z 

14a ; 

122 REM UNWRAHDLUNG IH EINE ALFHAHUNMERITSCHE "RRIRELE 

1E0 ZF=ITRErZ> 

17a ; 
134 G= 
124 FÜR I=1 Tü LENLZEI 
zag Dit Y/ALCNMIDELZE,-I.1>) 

21a NE=T I 

z28 3 

220 PREINT'QUERSUMME YÜOnN':Zz>"'=";0 
240 END 

FERD'T. 

LER ZUÜNrNNMIE 


EINGABE DER ZAHL? 1234598733 
QUERSUMME YOH 1223458793 =4 
rn m 
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2 EWIGER KALENDER 


Ein bekannter Algorithmus zur Berechnung des Wochentags für 
ein gegebenes Datum ist das Verfahren von Zeller. 
Für den Tag T, den Monat M und das Jahr J ergibt sich der Wo- 


chentag W aus der Formel 
W= [2.6M-9.2] + J + [J/4] + [H/4] - 2H + T 


dabei ist H die Nummer des Jahres innerhalb des Jahrhunderts. 
Die eckigen Klammern deuten auf die Anwendung der INT-Funktion 
hin. Zu beachten ist, daß Monate und Jahre nach dem spätrömi- 


schen Kalender berechnet werden müssen: 


Januar und Februar sind die Monate 11 und 12 des Vorjahres; 


März ist der Monat 1 und so fort. 


Der Wochentag ist natürlich nur mod 7 bestimmt; es bedeutet 


1) Sonntag 

1 Montag 

2 Dienstag 
usw. 


In Pascal können die Variablen Tag, Monat und Jahr zur Kon- 


trolle der Eingabe als Unterbereichstypen definiert werden: 
tag = 1..31; monat = 1..12; jahr = 1582..2999. 


Ein Jahr vor 1582 einzugeben ist nicht sinnvoll, da der Gre- 
gorianische Kalender damals noch nicht gültig war. 

Die Zellersche Formel wird als Funktion Wochentag definiert. 
Mittels der Mehrfach-Alternativanweisung CASE .. OF wird der 


gesuchte Wochentag ausgedruckt. 


In BASIC müssen die MOD- und DIV-Funktion durch die INT-Funk- 
tion simuliert werden, wie es beim Programm Quersumme gezeigt 
wurde. Zur Ausgabe des Wochentags wird am besten die Mehrfach- 
Sprunganweisung ON .. GOTO benützt. Damit keine Zeile mehr- 
fach durchlaufen wird, muß jede Zeile mit einer weiteren 
Sprunganweisung enden. Am Fehlen der GOTO-Anweisungen zeigt 
sich die Strukturiertheit der Programmiersprache Pascal be-: 


sonders gut. 
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ja 
ıi1u 
124 
128 
146 
154 
120 
170 
130 
136 
za 
21a 
22a 
220 
240 
PT 
Pad SIE 
era 
238 
2>a 
Su 
318 
3208 
Fa 
248 
Sa 
2E0 
ara 
Rt) 
320 
+00 
410 
420 
+30 
440 
470 
420 
470 


co 
= 


Wi in La BO BO 
we 
EEE 


Ewiger Kalender 


Praoarasm ewiger_ kalender input output), 

‘$ Jiezses Programm kerechmnet zu jedem Datum 
des gregar. Kalenders den WMochentzy 

mach dem Yerfahiren des Geirtli.Zeller#) 

t.pe  .„ahr 1222. .26033; 


monat = 1.„12: 
+24 = 1..21> 
IN BT) Jhr : Jahr; 
mar : meomat: 
+4 ı ta: 
c# %, 


function vachentasrtstagsmemanat?ı:sjahnrd sınteger ; 
var jr. ihdt:inteser; 
kein 
if m>Z2 then mı=m-2 
elze kein 
ms=m+1ä; 
.Jeı=3-1 
er: 
Jrs=) mad 104; 
Jhdt:= j div 100; 
wachentags=i.13#m-1>div Stijr div drjhldt Jdım d+air+t-ZHjhelt> mod 7 
end? 
# +) 
keyin C#Hauptprcaramm #%) 
writelnc "Gib Ta9g,Manat und ITahr ein!”>; 
readttg,.man,sjhr>e: 
write "der ",t9:2,"." ,maons2,.”.”,.J3hrsd.” ist ein ">: 
case wochentageta,meon,shr> cf 
Asurite ind "Sonntag" >: 
l:write ind "Montas”>; 
2swrite in: Dienstag”); 
3suritelnd Mittwoch” >; 
Jswrite inc" Donnerstag” >: 
Sswuritelnc"Freitag”>; 
Esuwrite Inc "Samstag" > 
end ı# case %#) 
erid. 


REM EWMIGER KRLENDER 
PRINTCHRFU1472 


REM HACH DEM VERFAHREN YOH ZELLER HIRD FUER JEDES DATUM 


REM DES GREGORIRHISCHEN KALENDERSE DER WOCHENTAG BERECHNET 
PRIHT" EWIGER KALENDER " 

PRIHT 

IHPUT" TRGS, MÜHRT , JAHR" ST.M,J 

IF 121532 THEH PRINT" JAHR VIERSTELLIG EINGEBEN" :SGOTO 1506 
PRINHT"DER";T3".":M>".":I2" IST EIN "; 


REM ZELLERSCHE FORMEL 

IF M>Z2 THEH M=N-2:6G0T0 258 

M=M+1G0:1=1-1 

H=INT« 1/1005 

J=T-100%&H 
M=IHTETF4SHINTEHF4D HINTER 1 3ENM- 1999 + Tr I-2EH 
IF za THEH WM=ul-7#INTCHT)SGOTO Zac 

IF Ya THEN YHeh+r:GOTl 236 


FRIMT"SZONNTRG":GOTO Jacı 
PRIHT"MOHTRG" :GOTO 4a 
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3424 PRIHNT"DIENSTAG"SGOTO 4a 
352 PRINT"MITTMOCH" :GOTO 4ası 
2er PRINT" DONNERSTAG" SGCTO 4a 
37&U PRIHNT"FREITAG"SGOTO das 
320 PRINT"SAMSTRAGS" 

230: 

+aa EHE: 
FKERD'T. 


EOIAJIGERrR HRALEHDER 


DER 24 .„ 12.1382 IST EIN FREITAG 


3 SIMULATION DER TEILERFREMDHEIT 


Das Nachvollziehen eines Zufallsexperiments mit Hilfe von Zu- 
fallszahlen nennt man Monte-Carlo-Simulation. 

Im folgenden soll die Wahrscheinlichkeit, daß zwei natürliche 
Zahlen teilerfremd sind, simuliert werden. 

Dazu werden mittels eines Zufallszahlen-Generators 19898 Paare 
von dreistelligen Zufallszahlen bestimmt. Durch Anwendung des 
Euklidschen Algorithmus wird der größte gemeinsame Teiler 
(ggT) ermittelt. Ist der ggT der beiden Zufallszahlen 1, so 
wird der Zähler der absoluten Häufigkeit weitergezählt. 

Die relative Häufigkeit stellt den Monte-Carlo-Schätzwert für 
die gesuchte Wahrscheinlichkeit dar. Man kann beweisen, daß 


hier die Wahrscheinlichkeit 


>. 8.6079... 
T 


ist. 
Da im Standard-Pascal kein Zufallszahlen-Generator definiert 
ist, muß ein solcher definiert werden. Dafür gibt es zahlrei- 


che Möglichkeiten; hier wird die Folge 


tx 1m+x 
n 


X = e - TRUNC (e n 


) 


mit einer zufällig gewählten Startzahl X, benützt. Sie ist im 
folgenden Programm als Funktion RANDOM gegeben. 
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Der Euklidsche Algorithmus kann durch eine Wiederholungsan- 


weisung formuliert werden: 


repeat 
rest :=m mod n; 
m:=n; n := rest 
until rest = @; 
ggT := m; 


In BASIC ist ein Zufallszahlen-Generator implementiert, der 


meist durch RND (X) aufgerufen wird. Dagegen gibt es keine 


Funktionsprozeduren. Der Euklidsche Algorithmus muß somit als 


Unterprogramm eingegeben werden. 





randamy=“ 
erde 


Kein 
repasat: 


rezts= m mad 
nıen:ne=re=t 
umtil rezst=d: 


garten 


a 
a 3er] 


= 


= 


ru 
* 


Kr 
u 


en 


e4in .#Haurteroar aumm 
zaehlers=ir (# Set 
=, Sreldide-al: Fr 
for ijs=l to anzahl 
kesin 
asetruncoktause 


vie 
. 


De u Er} 


De | 


Ksetrunee taus 
if kur za then 
gie 
haeufiakı=zaeh ler 
uritelne "ODie rel.H 
erid, 


RRERrLRRRLLL GGG LE LI a a Pi Pi fa Tu Pi Pi im m m m a 


Jia EU Ne DE NS) TON Bu e| ic 


in 


ai proaram teil lertremdhe 

1a !# PMonte-lzar lo-Simul 

zu da=zz 2 matuer liche Zahlen 

2a conzt anzahl = 

4 tausend = 

bl ri = 

15 BR FT 2o zaehler.i. a. 

va haeutiak:. s 

201 # 

33 funıtıion randemiırezal; 

ar Kaın 

1a KESSEL 

Pe z3ex-trum.iude 

3 

4 

za 
tunation satln.nsinteser) 
„a rezstrinteser; 


ıtiautput: 


” 
r 


Hahrzcheinlichkeit 
teillertremd zind #) 


aticon der 


} 


1a: 

1a: 

2, 193109268322 3 
integer ; 
real: 


sinteger: 


ri? 


* 
zen des 
> Start 
In {m 


er) 
Random #) 


-aehler 
Bei ae 


ndsr andom> 


end&r anıdam? : 


if watla,k>=1 then zaeh ler ı=ezachler+i 
“anzahl: 
zeufigakeit =". haeufiakıl a) 
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aa REM MÖHTE-CARLO-SIMULRTION DER MAHRSCHEINLICHKEIT 
11a REM DRSS ZWEI HRTUERLICHE ZAHLEN TEILERFREMD SIND 
120 PRINT"SIMULATIOHN DER TEILERFREMDHEIT" 
120 : 

14a Z=U0 

130 REM ERZEUGEH DREISTELLIGER ZUFALLSZAHLENH 
led FÜR I=1 TO 1a 

170 R=INHTEIBRAERNDOE1>) 

1220 B=INTEIARMIERHDIL2H>) 

139 GOSUB 2308 

zaa IF G6T=1 THEN Z=eZ+l 

ia NE=T I 

zz 3 

220 REM AUSGABE 

240 H=2.’1a00 

252 PRINT"DIE REL.HRELFIGKEIT=";H 

zei END 

era: 

224 REM UHTERPRÜGRAMM FUER G6GT 

230 REM EUKLIDSCHER ALGORITHMUS 

zag3 R=A-INTIR’E>HB 

IF R=&a THEH GGT=B:FRETURHN 

R=B :B=R 

2a GOTO >48 

FERDT,. 


= Im_L_HMT IcH DEF TEILEPFREMOHEIT 


„ 
m 
br! 


ul DO 
IKT LER 


[3 


AERARAE 


DIE FEL.HREUFIGEEIT= „ei: 


Lu 


4 CRAMERSCHE REGEL 


Mit Hilfe der Cramerschen Regel lassen sich lineare Glei- 


chungssysteme mit drei Unbekannten mittels Determinanten be- 


rechnen. 
Beispiel: _ _ 
6x, AxX, + 3X, = 7 
3x, + 3X, - 8X, = -1]1 
2x, + AX, = 14 
Die Hauptdeterminante D_, 
-4 3 
3 5 -8 = 202 
2 ) 4 


verschwindet nicht und zeigt so, daß das obige Gleichungs- 
system eindeutig lösbar ist. Setzt man die rechte Seite des 
Systems für einen der Spaltenvektoren ein, so erhält man drei 
weitere Determinanten, aus deren Wert sich die Unbekannten be- 


rechnen lassen: 
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7 -4 3 
_ 1 |_ _g|_ 222 _ 
x = D_ 11 5 8| = 302” 1 
14 6 4 
6 7 3 
I _ _g| - 494 _ 
y- D_ 3 11 8| = 307 7 2 
14 4 
-4 7 
_ _ı _ _ 696 _ 
2=5 5 11| = 302 3 
0 2 e) 4 


In Pascal läßt sich das Einsetzen der Spaltenvektoren durch 


Wahl eines passenden Variablentyps direkt programmieren: 
type vektor = array [1..3] of real; 


Auch zur Determinanten-Berechnung kann eine geeignete Funktion 
definiert werden. Damit das Feld der Koeffizienten nicht ko- 
piert werden muß, wird das Feld - wie auch in den folgenden 
Programmen - als globale Variable behandelt. 

Da BASIC keine Funktionen mehrerer Variablen kennt, muß die 
Determinanten-Berechnung in einem Unterprogramm vorgenommen 
werden. Dazu muß allerdings die Variablenübergabe explizit 
programmiert werden, da BASIC nur über globale Variablen ver- 


fügt. 


1A program cramersche_rewse \le.coutput): 


ia canzt n=3; .# Zahl der Unkekannterns#)> 
120 tope vekter=zarrzu[.1l..n) of real; 
122 „ar a,.k,0,d : weht: 
14 hauptdet,.=z,v,.,.zı real; 
indes : 1..rnr 
Die 2 %) 
+uncotıon Jdeterminmantern ,‚ z,twwelktordıreal, 


kein 

determinentes=er [1]#=2[2J#t031J+rC2J#2:[2J#tCo1I+rl3lsesliletle) 
-" [3J#=2[2J#t[1I-r[23#=[11JI#t023I-r [1JI#E[3 Ir LEI 

erd:c$ determinante %#%) 


Deo #) 
Pracedure eingake: 

kesırı 
all):sE,.A;k[E1l)se-4.M;cl1]Jıe2.8:0[01):=t. 0; 
alz2J,=23,.2:k[2JI,=i.B;c[2),2=-2. Brdl2]se-11.0: 
al?)Js=2,. u; kl) sei.ü;:cl2)eed.Brdl3):eid. 0 


er? 

CH #5 
keuın © Haupteroaramm #. 
ejinmake 
Kauptdet:=deterninaentela,k.: 
erite ine" Hauptdeterminante- 
it =@kzs{ihauptdet)/1l. Bein 


BR 


Lu u 


hauptdget); 


WdE PERS BR N NS Di Mer m nm m 


Eifer gig ini Ba ic in 
DE a u Sc or cv a 5c Bor Br BR 9 BB Sc BR SCH BE 39 1 59 Sc 1X 9 Sc 19) 
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then writein. "Gleichungszmatrix zingulser”) 

else keyirn 
zs=determinanteid,k,c’chauetdet? 
„s=determinantela,d,c’rhauptdet; 
z:=determinenter a,b. d’rhauptdet; 
writelni"z=",2," va7.y,” Ze’, 


erc 


SE 
GAn5mEh 


De a Su Er Ey Eee EEE 


end, 





FEN CRAMERSCHE REGEL FUÜER 3 UNBERRHHTE 
FEINTCHREN 1472 
PRINT"CERANERSCHE REGEL" 


e, 
[2% 


DIM Ac3,45,.Di3,.45 

FOR I=1 TO 3 

FREAD ACI,1D.RCI-2D.AEI,ED.AUL.4) 
HEST I 

REM HAUPTDETERNIHAHTE 

FÜR I=1 TO 3 

FÜR J=1 TO 2 

: DEI.DD=ACI,D 

HEST I 

HEST I 

GOEUB 54a, 

Da=p 

IF Da=a THEN FRINT"GLEICHUNGSETSTEM NICHT EINDEUTIG LOESBAR" SEND 
REM #-DETERMIHANTE 

FOR I=1 TO 3:D@1,.1>=ALI,.45oNEHT I 
GOSUB 54a 

Def) 


FEN T-DETERMINAHTE 

Für I=1 Tü 3 

; Dr. I,.15=HC1.15:DE I.Z>=nn 1.dn 
HEAT I 

GOSUB Sg 

D'yr=Q 

REN Z-DETERMINAHTE 

FOR I=1 TO 2 

: DE I.Z>=NeEl1,23:DE1.2>=RK 1.4) 
HEA=T I 


GES SIE EZ 


DD LOS CHEN Bi Re N Bin Bu Ri ib Min Bi Ar a I a IB EN a Pi 


DIN IRB BR BR BL BB RL LS SS N EI DI DM DB MD N TDI N DM N Fe rennen nn 


a GÜSLB 54a 
a DEz=D 
[ts Bu 
a REM AUSGREE 
a PRIHT"LOESUNG" 
a FRINT"S=" O8," Te" ODE," Ze DD 
ia END 
Ps 
2a BEN UNTERPRÜÖGEAMM FUER DETERMINAHNTENBERECHHUNG 
353 OD=071,.1>#D0%.2,2:#D0.2,3>+0D7 1,2530. 2.2>%4003,.15+D7C1. 324002, 13$003.20 


nn 


za D=D-Di1,33530°2,2>3D052,195-DE1,2>#%0C2,1>#00.3,3>-DE1, 15002, 29ED E25 
seRn RETUFMN 
bu 5a u 
sag DATA &,.-4,.2,7 
sag DATA 2,.5,-2,-1l 
san DATA Z,4a,.4,14 
PRERAD'T. [FRPMNERSCHE REGEL_ 


a3} 


.! 


LOESUNG 
“=1 Y=2 


Mi 
u 
1 
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5 KOMPLEXES RECHNEN 


Die für viele Anwendungen benötigten komplexen Zahlen lassen 


sich in Pascal bequem als Verbund (RECORD) vereinbaren: 
type komplex = record re, im : real end; 


dabei stellen re und im den Real- bzw. Imaginärteil der kom- 


plexen Zahl dar. 


Die Addition und Subtraktion der komplexen Zahlen ist kompo- 


nentenweise erklärt: 
(atbi) # (c+tdi) = (atc) + (b+td)i 
Die Multiplikation ist definiert durch 
(atbi) (ctdi) = (ac-bd) + (ad+tbc)i. 


Die Division kann als Multiplikation mit dem Kehrwert 


(c+di)”' = ——_ (c-di) 


c?+a° 
durchgeführt werden, dabei stellt der Nenner den Betrag der 
komplexen Zahl dar. 
Alle Rechenoperationen können in Pascal als Prozedur erklärt 
werden; die Ergebnisse werden als Variablen-Parameter ans 
Hauptprogramm übergeben. Die Betragsrechnung kann als Funktion 
definiert werden, da hier ein reeller Wert erhalten wird. 


In BASIC gibt es den Datentyp RECORD nicht. Die komplexen Zah- 
len werden deshalb als reelles Zahlenpaar aufgefaßt. Wie in 

den vorausgegangenen Programmen könnte man die Rechen-Prozedu- 
ren wieder als Unterprogramme schreiben. Dies bringt hier aber 
keine Vereinfachung, da diese Programmteile nur einmal durch- 
laufen werden. Dagegen wird die mehrmals benötigte Ausgabe als 


Unterprogramm formuliert. 
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13 pragrzaun komplexes _ rechnen output>?; 

119 tupe komplex = recard 

124 re.imsreal end; 

124 var a,kb,.cı Komplex. 

i4a rocedune zuzsgakerzckamg lei)? 

159 kesin 

lied it x. im then writerz.re,”+ti#" 223,2. im) 

1” elze writelz,re,"-i#" 22, ahzix. im>>} 
wmitein 


fe 
E: 


aa 
132 endsi# at zusgakex) 
za prucedure zummecz,vskomp lex zvar zıkomp lex); 
zia kein 
Pan zZ, ressx,retwv.r©e; 
za z. ims=x. imtw. im: 
24a en; 
242 procedure diftferenztix,v komp lexsvar zikomp lex>; 
244 besir 
z4E zZ. reszsx,re-v.re; 
” z. ims=x. immw. im 


43 erg 
ad ErocLedure praoduktiıx,vosköomp lex zvar zikomp lex) 
keyir 
Zz,res=2,refv,re-x. im. im? 
z. imssx.re#t,. imtx. imfv.re: 
ru} 
fur-tion betraylxskomp lex?’ sreal: 
keusir 
ketraygszseartisarz.rebtzarlz.im>) 


er: 


an Bu 1 BE Ua BE U ED 
Nu En Eu 
ET er a ce Bu A Eu oe u 


x) m 


Sal Ka da 


Echt pa 1) 


DE 


249 peraLedure inmwerzez. x :konp len mvar zıkomp lem); 
Ga begin 
ZEN z,.reı=sx,resarketraulx)>} 
ar z. ims=-x. im’ za betragt s)) 
3209 erd} 
N erocedure eirdake 


Agı 
u EN LE | 
Do) 

dt 

® 

+ 

a“ 


kesirn 
e;= rn MET 

Fa K,re:=5,0:k,ims=2,.0 
420 erıdz 

444 beyin <#Hauptproaramm #5 
445 ejrgake + 
420 zummera,.Kb,.0d} 
4 wumiter "Zumme="> 
472 Jifferenzia Ben 
4,7 umjiten "Ditter nz 
va Produktia,kb,.che 
uritel "Produkt=""rausgakbelc)he 
irversesik,c); 
HERR 
uriter "Quotıient=") ausgabe: 


Da Ba ı Bu et Su 


Pie a ich ı 
DC ct Bu ca Er Br Ey 


1 2 [De | a 


mm nn en 
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- 
Id 
Es 


FEN A GRUHDRECHEHRRTEH FUER KÖNMPLE=-E ZAHLEN 
PRIHNTCHRE14TD 

FRIHT" K.ÖMPLE=ES FECHHEH"SFRIHT 

DIMmM REiUZ2>.IMGZ> 


ER) 


FEM EINGABE 
FERD REr1:,.IME1>,.REr2.. IM 


So 


DO 1! 


REM ADOITIOH 
RE=REÜL>+RELDD 
IM=IMi1>+Ir.2> 
PFEIHT "SUNmE=": 
GOZUE 452 


En 


FEN SUBTRAHTIOHN 
FE=RECI>-FENZED 
ImM=INe1>-ImMdz" 

FREIHT "DIFFEREHNZ=": 

GOSUEB 450 

FEN MULTIFLIERTIOH 
RE=FEr:.I>#FErLZ2>-IMG15HIMG 25 
IN=RELI>HIMCZ>H+REL Z>HING1D 
PEIHT "PRODUET="; 

GOSUE 458 

REN DIWISIOH 
EETRRG=SRELE>TZ+INGZDTZ 
FE=FEÜLD#RELZIHINGEIDHINGCED 
IM=-FEe1>#IMLZ>+REL2>#IME 15 
RE=FE.-'BETRARG : IM=IM.-EETRRG 
410 FRIHT "OUOTIENHT=";: 

42a GÖOSUE 450 

330 EHD 

44: 

450 REM UHTERPRÖGRANMM ZUR RUSGREE 
350 IF Imza THEH FRIHNHT RE:"-";RESCINMD "HI" SGOTO 4a 
47a PRIHT RE:'"+":IM> "HI" 

430 PRINTSRETURH 


Ja il Ka Li Pa Pa Pi TE P5 To Pd Pi Pa m Dt 


Zn hin Bi MR IN Bu en a ih aD id Da 


a a 


42a 5 
sagt DATA 7.-2.3,-3 
RERD'T. 


H-OÖOMPr_ En E> FECHHNEH 


SUMME=S12 + 1%1 


DIFFEREHZ= O2 O-I#1 


PRÜDUKT= 41 +11%I 


QLIOTIENT=. 352941176 - . D11rERcHl 





6 PRIMZAHLSIEB DES ERATOSTHENES 


Ein bekanntes Verfahren zur Bestimmung von Primzahlen ist das 
Primzahlsieb des Eratosthenes (ca. 280-200 v. Chr.). 
In der Menge {2,3,4,5,.....,N} werden zunächst alle Vielfachen 


von 2 - jedoch die 2 selbst nicht - gestrichen, wie man sagt, 
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"ausgesiebt". Sodann werden alle Vielfachen der nächst größe- 
ren, im Sieb verbleibenden Primzahl entfernt - jedoch :nicht 
die Primzahl selbst. Das Verfahren setzt sich in der angegebe- 
nen Weise fort, bis die Vielfachen einer Primzahl gestrichen 
werden sollen, die größer als VN ist. Die im Sieb verbliebenen 


zahlen sind genau die Primzahlen der gegebenen Menge. 


In Pascal kann das Sieb als Mengenvariable vereinbart werden: 


sieb : set of 2..n; 
Folgende Mengenoperationen sind definiert: 


A+ DB Vereinigungsmenge 
A* B Schnittmenge 
A- B Restmenge, 


ebenso folgende Booleschen Operatoren 


IN Prüfung auf Element-Eigenschaft 


<= Prüfung auf Teilmengen-Eigenschaft. 


Die Abfrage, ob ein Element der Menge angehört, kann somit 


formuliert werden als: 


for zahl := 2 ton do if zahl in sieb then ... 
Entsprechend kann das Aussieben durch Bilden der Restmenge 
sieb := sieb - [vielfach*zahl] 
durchgeführt werden. 


Genausowenig wie den Typ RECORD kennt BASIC den Datentyp Menge. 
Man kann sich aber wieder mit einem Feld A(I) behelfen; gehört 


die Zahl i zur Menge, so wird A(I) = 1 gesetzt, andernfalls 


zu Null. 
Das Sieb wird gefüllt durch 
FOR I=2 TON : A(I)=1 : NEXT I. 
Eine Zahl j wird ausgesiebt durch 
AU) = 2. 
Die Ausgabe der im Sieb verbleibenden geschieht durch 


FOR I=2 TO N : IF A(I)=1 THEN PRINT I;: NEXT I. 
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Wie man sieht, ist die Programmierung mit Mengen wesentlich 
eleganter. Ein Nachteil des Mengentyps in Pascal ist jedoch, 
daß Mengen nicht ausgedruckt werden können. Dasselbe gilt 
auch für Unterbereichstypen (siehe Programm Ewiger Kalender) 
und Verbunde (siehe Programm Komplexes Rechnen). 

Ein zweiter Nachteil des Mengentyps ist, daß die meisten Com- 
piler die maximale Elementeanzahl stark einschränken. 


Der hier verwendete Compiler läßt nur 127 Elemente zu. 





lag prouram eratcsthenestautput>: 
lila ı# Frimzahlzsiek des Eratosthenes #5 


120 const n = 1297 
130 var ziek zet of 2,.n: 
140 i-zahl,vielfach s integer; 
120 CH 4) 
lEO procedure zuzsyakı: 

IT a) kiinteyer; 

keyirn 


for kı=z ton ca 
if kın ziek then writeik:d) 
ends# zuzaake #) 
De © #) 
kein # Hauptproaramm #) 
siek;s= [2..n]; 
for zahl:=e2 to truncisartin>s> do 


kesirn 
if zahl in siek then 
kbeyin WFzauszsieben#) 


„jeltrachs=z; 

while vieltrach {= nn div zahl de 
kegyın 
zjeks=erjek-[vielfachtzahl]s: 
„jeltfachsev3jelfach+ti 

end 


zz a a a a a a a 


VUN Ta N Ta fa N 
De Dig Thin Bis Bd 
DR) 


Ki ya Kai El Ka Ka Kar d 
De BE Du BR 
IKT) 
T 
I 
zn 


a 
ni 


KT uH 


a z 


i 


er. 


jan REM FRIMZAHLEIER DES ERRTOSTHENES 
11a FRINTEHRE 14T 

122 FREINT"SIEE DES ERATÜSTHENES" SPRINT 
1:0 3 

i4G IHFUT"OBERE GRENZE" HN 

152 oIM FiH> 

12H ; 

1’ıa REM FUELLEN DES SIEEN 

123 FüR I=2 TO MH 

134 ; Fil>b=1 

zaa NEAT I 

li: 

zz REM AUSSIEBEN 

20 Für I=z2 Tl ZUR.ND 

za 5 IF Fiıld=a THEH 250 

zsA s T=1#] 

zei: FÜR K=T Tl H STEP I 

zvas : Pike 
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zo 2 HMEAT K 

z30 NE=T I 

B.1 515 Bu 

10 REM AUSGABE 

z2a PRINT"FRIMNZAHLEN EIS" :Ne"SIND" 
224 Für I=2 TO MH 

2409» IF Fel>=1 THENH FRIMT I: 

=aa NEAT I 

250 END 
FERD'T. 


Mr ImzZrAHL_=>IEF OE:2 EFRTOSTHEHEDZ 


ÜBERFE GRENZE ? 125 


DIE PRIMZAHLEN EIS 125 SIND 
ri 
1032 1: 


a If 


1: Pa 
i 
i 

IT —- 


on 


ca 2 
Deu eg = 

Zi 
Il a m 
u 1 u! 
ka ii 
Pan 

[N m 
u Et 
u a 8 

-j dd 
et u 4 
Li Li 


ib 5 
EN N 
I Ny m 





7 BINÄRES SUCHEN 


Als Beispiel eines rekursiven Verfahrens sei das binäre Su- 
chen in einer geordneten Liste behandelt. 

Von der gegebenen Liste mit bekannten Index-Grenzen wird zu- 
nächst der mittlere Index bestimmt. Vergleicht man das zuge- 
hörige mittlere mit dem gesuchten Element, so wird entweder 
das gesuchte gefunden oder aber diejenige Hälfte der Liste 
bestimmt, in der es sich befindet. Der Suchvorgang beginnt in 
dieser Hälfte von neuen. 

Das Verfahren endet sicher, da die Anzahl der zu durchsuchen- 
den Elemente jeweils halbiert wird. Spätestens beim letzten 


Element entscheidet es sich, ob die Suche erfolgreich war oder 


nicht. 


In Pascal läßt sich dieses rekursive Verfahren direkt program- 
mieren: Die Prozedur binärsuche ruft sich bei jedem Halbie- 
rungsschritt mit veränderten Grenzen selbst auf. Die Variab- 
len-Parameter index und ergebnis übermitteln das Ergebnis der 
Binärsuche. Die Suche endet spätestens, wenn die zu durchsu- 
chende Liste nur noch 1 Element enthält, da bei einem erneu- 
ten Prozeduraufruf der obere Listenindex kleiner als der unte- 


re wird. 
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Da es in BASIC keine rekursiven Prozeduren gibt, muß das binä- 


re Suchen iterativ programmiert werden. 


Bezeichnet I,J,K den Index des kleinsten, größten bzw. mittle- 


ren Listenelements, so wird I=J gesetzt, wenn das gesuchte 


Element sich in der linken Listenhälfte befindet, andernfalls 
K=J. 
Das Programm endet spätestens, wenn die Liste ein Element ent- 
hält; d.h. wenn K-I=]1 gilt. 





ac 
1a 
124 
132 
140 
120 
120 


J 


a 


Kr 
a 


“ = Fa 
Zn 


TI a N DE De Pa a 
1 


LRIRERLH IR IB IN IT EB BR RB BR LBS Ga LG id LE BO RE POS Pi RI PS DI RE m 


OU NT a ne N BE u 


_ 
ih 


Praaması Kinzerzuchetlautput>; 
camz=t HM = 1; 
tepe lızste= arraurlid,.n)] or integer; 
war zahl,.ztelles integer: 
gefunden : kaclearız 
1 : liste; 
# + 
Pracledyure Kinszerzsuchelunten,.caken, zsinteger zoar indexsinteger: 
„ar ergekmızıkaoclean.? 


„var mwittesinteger; 

kein 

nitter=eickentuntend div 2: 
ıt cken % unten then Kesırn 


ergekmnizı=stfalse; 
index ;=u 
er 
el2e 
ı* lEmittelsz then 
keyin 


ergeknizs=true?r 
index ,=mitte 
ericd 
l= 
ı* lEmittel>z then 


ıYı 
ıYı 


Kinszserzucherunten.nıtrte—-il .,. zahl. =telle,. gefunden. 


elzE 
Kinserzuchernitte+rl . cken. zahl. stelle. se funden. 
end: !* Kimnserzuche #5 


3 *. 
pPröücedure eirgahe; 
keairn 


zn tierte Liste %) 

1[1]:=12: 1[L211s=23H+ 1[L2]J 2827: 
1[43s:=32: 1[5312:=35; 1[e}]smde; 
1[’]s= 9, 1[2]ı=06: 
1[.31:=31:1[.120]:=33 

end ı#$ einsake #5 


kesin ©$ Haupterouramm #95 
eimsahe: 

writelnd "gikb gesucht 
readlzahl>; 
Kinaerzuchei.n,zeh |l,ztelle, getunder>?: 
ıf aefunden then 


ıl ein !"D: 


ıyı 
Ts 
u 


0 


wmitelnmi " Gezsuchtee Element an” ,.=tell 


lz 


ıYı 
ıı 


.. 


ı 


Zr stelle” 


wıte lm "Gezuchtes Element mıcht setunden” > 


rc, 
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iaa REM BINRERES SUCHEN 

ila PRINTCHREEIYHTD 

2a PRINT " BIHNRERES SUCHEN" s:PRIHT 
120 : 

14G REM EINLESEN 

1229 RERD N:DIM ACH» 

FÜR I=i TÜüH 

RERD ACID 

HE=T I 

IHFUT"GESUCHTES ELEMENT" 28 3UCH 
PRIHT 

REM AHFRHNGSMERTE DER INDEAS-GREHNZEH 
I=M:K=H+1 

REM SUCHSCHRITT 

J=IHT%CI+K)2> 

IF »5SUEH=ACT> THEN PRINT"GESUCHTES ELEMEHT RAH":JT:". STELLE" sEHD 
IF »SUCH > REIT THEH I=I:G0T0 31a 
=.T 

IF K-I=1 THEHN FRINT"GESUCHTES ELEMENT NICHT GEFUNDEN" END 
GOTO 272 

EHO 

REM DATEN 


DATA 12,.22,27,22,35,94,090,686,31.33 


M io do nd di BG RI NED NS RE DE DI DI DE m m m 


Leo in Bis Re DIT in Bi Bd 
OS Sun N US GH GE EEE aa 


El 


EIHRERES SUCHEH 


GEZSUCHTE> ELENEHT? ©6 
GESLUCHTES ELEMEHT AH 2.STELLE 


8 INTERVALLSCHACHTELUNG FÜR NULLSTELLEN 


Eine wichtige Anwendung des binären Suchens ergibt sich bei 
der Nullstellenbestimmung durch Intervallhalbierung. 

Nach Eingabe eines Intervalls wird durch Berechnung des Funk- 
tionswerts an der Intervallmitte entschieden, ob sich die ge- 
suchte Nullstelle in der linken oder rechten Intervallhälfte 
befindet. Sodann wird das Intervall auf die entsprechende 
Hälfte reduziert. Setzt man das Verfahren in der angegebenen 
Weise fort, so wird die Nullstelle durch immer enger werdende 
Grenzen eingeschlossen. Um sicher zu gehen, daß die Funktion 
im eingegebenen Intervall eine Nullstelle besitzt, wird die 
Funktion auf Vorzeichenwechsel geprüft. Denn nach dem Zwischen- 
wertsatz für stetige Funktionen hat jede stetige Funktion in 


einem Intervall, in dem sie das Vorzeichen ändert, mindestens 
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eine Nullstelle. Obwohl danach die Nullstelle sicher gefunden 


wird, könnte es sein, daß die eingegebene Genauigkeit nicht 


erreicht wird. 


ausgegeben. 


In diesem Fall wird eine entsprechende Meldung 


Analog wie bei der Binärsuche wird das Intervallhalbierungs- 


verfahren als Prozedur intervallhalb definiert; sie ist jedoch 


iterativ formuliert, da iterative Verfahren meist schneller 


als rekursive sind. Die Intervallhalbierung endet, wenn die 


gesuchte Genauigkeit erreicht worden ist oder mehr als 35 Ite- 


rationsschritte benötigt worden sind. Die Variablen-Parameter 


mitte und ergebnis übergeben die Ergebnisse der Nullstellensu- 


che an das Hauptprogramm. 


Das BASIC-Programm entspricht wegen der Nicht-Rekursivität dem 


Pascal-Programm; jedoch sind alle Variablen global. 


nen ne 


1a 
ı1a 
128 
134 
140 
150 
128 
ira 
130 
1>8 
eu 
Pr 
zei 
ea 
40 
ze 


ei 


.e 
'e 


} 


1 


Je la di WE GE LEI LE LE a BE Pi Ra Ti 
ee 5 a a a ee NE ET er 


BDc® ce cr Ba ce cn BE cn Er EEE En 


L 


J 
- 5 
m 5 


420 
420 
440 
aa 


Fraouaram ıntervallhalkierunslinput,.output>: 


ar nullstelle,z,'w ‚real; 
varzeicherwechzel,getunder cskaoc learn; 

u #5 

function funktıonixerealbsreal: 

kesir 

funkticon:sssarcx)-2.0 

erde 

Ei # 

eracedure intervallhalkl function Fıreal:a,kır 
Te 

conzt epzi laon= 1.Ge-KHE; 

ur zaehler: integer; 

heuirn 

zaeh ler :s=ü; 

repeat 
nittes=lark)/’2.0; 


zaeh ler s=zaehler+1: 
if Fla)s$tlonmitte)>in then aszsmıtte 
else kı=mitte 


until Caksck-a)T epsi laon>arczaeh ler.35, 
ergeknisı=lizaeh ler i=23) 3; 

end? 

Ci 4) 


kesin C#Hauptprcar zum #%) 
umiteilne "Gib Intervallgrernzen ein!”) 
readli,.w)3 


varzeicherwechzels=ctunktioncx>$Ftunktionlvbils 


if vorzeicherwechzel then 
ein 


ealroar mıttesresl: 
er gekmizskoolean)? 


intervallhalke funktion, x.v.mullstelle,getungden)?; 


ıf gefunden then writelnc"Hullstelle „ul late L leere 
else write ln "Genzuiskeit nıcht erreicht” 
ericd 
else writelnd "kein Yorzeichermechzel”>: 
er. 
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iada PEM IHTERVYALLHÄALRTERUNG FUER HULLSTELLEH 

1id PRIHTCHREK 1STD 

124 PRIHT" INTERY. HALBIERUHG FUER HULLSTELLEH'":PRIHT 
1234 3 

id4d REM EINGABE DEF FUHRTIOH 

17a DEF FH -Z 

lied > 

174 IHPUT"IHTERYALLGREHZEH RIE"+R,E 

120 EFS=1E-&:REN REL.FEHLER 

Sa 5 
IF FHFEIR>AFHFIE> <=a THEH 24a 
PFRIHT"EKEIH YVÜRZEICHEHMECHSEL " sEHO 
REM INTERV/ALLHALSTERUHG 
M={HA+B).’2:REM IHTERYALLMITTE 

IF FHFIR>SFHFIM> SG THEH BeM:sGüTo Zr 
A=hi 

IF ABEBSCB-A> > EPS#RBSCH> THEHN Zd4d 
PRIHT:PRINT"HULLSTELLE=" se A+BI 2 

EHO 

y, 


TEEN G 
Bagamaaanmm 


IM 10 Po Ro Po Po Po Po Po Pa Pl Ri m 


r 


27 


IH TER“. HA_EIERFRULHGS FUErR ÖAOHUÜlLLD=TEILLEH 


HALLSTELLE= 1.41421256 


9 ACHT-DAMEN-PROBLEM 


Ein Verfahren zur Lösung eines Problems, das durch Probieren 
vor sich geht und bei dem Schritte, die in eine "Sackgasse" 
führen, wieder rückgängig gemacht werden können, heißt Back- 
tracking-Verfahren. 

Ein solches Verfahren soll nun zur Lösung des bekannten Acht- 
Damen-Problems benützt werden: 8 Damen sollen so auf einem 
Schachbrett plaziert werden, daß sie sich nach den Regeln des 
Schachspiels nicht bedrohen. 

Wegen der Vielzahl der Lösungen ist das Problem kaum von Hand 
zu bewältigen. Es gibt nämlich 92 Lösungen, von denen jedoch 
nur 12 nicht durch Spiegelung oder Drehung auseinander hervor- 


gehen. 


Im folgenden Pascal-Programm wurden die Variablen wie folgt 


codiert: 
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Bedeutung 
x [il =75 Dame i steht in der Spalte j 
a [j] wahr keine Dame steht in Zeile j 
b [k] wahr keine Dame steht in der Hauptdia- 
gonalen oder einer dazu parallelen 
c [kl wahr keine Dame steht in der Nebendia- 


gonalen oder einer dazu parallelen 





Wie man der Skizze entnimmt, sind die Hauptdiagonale und ihre 
Parallelen durch gleiche Zeilen- und Spaltensummen gekennzeich- 
net; entsprechend die Nebendiagonale and ihre Parallelen durch 


gleiche Differenz aus Spalten- und Zeilenzahl. 


Setzen einer Dame kann so formuliert werden: 

x[i] := j; alj] := bli+j] := cli-j] := false. 
Entsprechend bedeutet 

alj] := bli+tj] := c[li-j] := true 
das Entfernen einer Dame vom Brett. Eine Dame kann gesetzt 
werden, wenn gilt 

aljJ and bl[li+tj] and cli-j] = true. 


Die Prozedur versuche setzt durch Probieren Damen aufs Brett: 
Gelingt dies ohne Bedrohung durch eine andere, so ruft sie 
sich selbst auf und setzt die nächste, andernfalls nimmt sie 
die Dame vom Brett. Sind 8 Damen plaziert, so wird die ent- 


sprechende Stellung ausgedruckt. 


53 


Acht-Damen-Problem 


In BASIC muß das Setzen wieder iterativ formuliert werden. 
Der Zeilenindex i gibt gleichzeitig die Anzahl der bereits ge- 
setzten Damen an. Weiterzählen von i bedeutet Setzen einer 
weiteren Dame, entsprechend wird i vermindert, wenn eine Dame 
vom Brett genommen wird. Die Bedrohung einer Dame wird wie 


folgt geprüft: 


FOR K=1 TO I-1 
IF D(I) = D(K) OR ABS(D(I)-D(K)) = I-K THEN ... 
NEXT K 


Das Weiterrücken einer Dame wird durch 
D(I)J) = D(I) +1 


erreicht. Sind 8 Damen plaziert, so erfolgt ebenfalls die 
Ausgabe der Damenstellungen in graphischer Form. 


nm 


120 eragram achtdamerprcok lemloutput>?; 
11a vor i:integer; 


124 sıarrawl[l1l,.2] of koclean;ü#keine Dame in Zeile if) 

1:20 ksarraw[l2..1&] 0 kocoleansl#keire Dame in Hauptodiag. i#) 
140 cıarrawrl-7..7) of kocoleanzc#keine Dame in Hekendiag. i%> 
154 ssarrarlil..2] © integer:l#Ffoasition in der Spalte i%#) 
1e0 CH #) 


1vG procedure zausgaker 
120 vor kıinteger?; 
130 begin 


za for kı=l t0 © go 

zie urmiteixz[k)Jsd4>r 

zz wmiteiln?z 

220 ende zuzsgakei) 

24a #> 
250 procedure verzucheri :ınteser >} 

Ze san jsinteger: 

2va kbesin 

zo For js=l tu © oda 

za it =[j]l and klı+j)] and clı-j) then 
za keyirn 

210 »Lils=i?: 

aa alj)ı=talzseskli+ti])sstalse: 

22a -[li-ijils=talse; 

24a it i«@ then verzucheritri> 

za else zusgakhe; 

ts) sLjJs=true:kli+rils=etrue; 

ara 2Li-jJ]Js=true 

a2 rd 

230 endzı# verzuche#) 

ad vH +) 
dig kesin !* Haupteraoaramm %) 

42a tor ı:=1 to 5 do alılsetruer 

422 +cor ı3=2& to 1E& do klil:sstrue:r 
Ts Fo ‚=-7” tu 7 da clils=etrue; 
+ verzgchenl‘ 

Jed er, 
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U 1111111 nn mn nn 


FE 2-EOHMEH-FRÖOBLEM 
FREIHNTCHREEIHTO 
bImM Di) 


Da Bat Sa a a RE u 
U NT 19 Se u 


A ;FEMN ZEILEHINDER 


= 
De 


1 
=] THEN zZ20 


- - 

pr che 

ne | rD nn Ti m 
N 
un BE u | 

u Be 


DE ET 


HE=T E 
IF 12 THEHN 1rE 


Pi Ri fa Dame ra 
06} 


1 Pl dh 


tz 


. 
Dr 


Für KE=1 TO 


Go lo 


1» REM HEUE DANE AUF 


FREIHTCHREE1I4TOGREM RUSGREE 
Z=Z+15;PREIHT Z>". LOESUNG" GSPRIHT 


REM ZREHLER FUER LORTUHNG 


BRETT 


FÜR k=1 TO I-1 :FEEM SPRALTEHINDER 
IF CXl»=DrcE> OF REBESCDKT>-Dcko.=I-h THEN 220 


= —- OCOHIEH - FRÖBL_ Er 


z2a FÜR L=1 Tl 1.LOEFSURG 
zaa IF Dikb=L THEH PRINT" ";:60T0 Zu 
aaa PRIHT"EE "5 

31a NE#T L 

224 PRIHNTSFRIHT 

322 HEAT K 

Bi 15 Bu 

>52 I=I-1 : REM DAME Yon ERETT 

aeg DriIo=oil:+l :REM DRANE FUECHT YOR 
aya IF DeIsl=2 THEHN 1984 

222 I=I-1 :FKEM DANE YO ERETT 

334 IF I>a THEN Zei 
REAL. 


10 LOGIK-AUFGABE 





Durch Rechnen mit Wahrheitswerten können Logikaufgaben wie 


folgende gelöst werden: 


3 Freunde A,B,C kommen unter folgenden Bedingungen zu einer 


Party: 


(1) A kommt nicht, wenn B nicht kommt 


(2) B und C kommen gemeinsam oder gar nicht 


(3) A kommt, wenn C nicht kommt und umgekehrt 


In Pascal können Wahrheitswerte (BOOLEAN) wie folgt verknüpft 


werden: 


A oder B durch 
und B 

impliziert B 
äquivalent B 


xor B 
ausschließ.Oder) 


> >» m » 


a 


a 


or b 
and b 


a<= b 


a 


a 


=b 


<> b 


Die obigen Bedingungen sind erfüllt, wenn 


(not b<=e not a) and (b 


c) and (a <> c) 
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wahr ist. Die entsprechenden Wahrheitswerte werden aus der 
Menge aller Kombinationen ausgewählt, die durch 3 verschach- 


telte Zählwiederholungs-Anweisungen von der Art 
for a := true downto false do 


erzeugt werden. 
In BASIC muß das Rechnen mit Wahrheitswerten durch arithmeti- 
sche Ausdrücke simuliert werden. Folgende Codierung wurde ge- 


wählt 


wahr entspricht +1 

falsch -1 

nicht A -A 

A und B sgn (A+B-1) 

A oder B sgn (A+B+1) 
A impliziert B wenn A<=B dann 1, sonst -1 
A äquivalent B sgqn (A*B) 

A xor B -sgn (A*B) 


Die Wahrheitswert-Verknüpfungen werden in getrennten Unter- 
programmen durchgeführt. Der obengenannte Lösungsterm wird 
schrittweise berechnet, indem das jeweilig benötigte Unter- 
programm angesprungen wird. Die dazu notwendige Parameterüber- 
gabe wird mit Hilfe der Variablen P und Q durchgeführt. Analog 
zum Pascal-Programm werden alle Kombinationen der Wahrheits- 


werte mittels dreier verschachtelter Schleifen 
FOR A=1 TO -1 STEP -2 usw. 


erzeugt. Diejenigen Kombinationen der Wahrheitswerte, die den 


Wert 1 ergeben, werden ausgedruckt. Hier ergibt sich 
A=-1 B=1 cC=1I; 
nur B und C kommen zur Party. 


In einigen BASIF-Tialekten - wie z.B. bei Commodore - sind 
die Booleschen Operatoren AND, OR und NOT definiert. 


-1 AD O = O0 


Es 14ßt sich zeigen, daß diese Verknüpfungen der Zahlen C und 
-1 den Gesetzen der Aussagenlogik entsprechen, wenn man 

-1 als "wahr" und C als "falsch" auffaßt. 
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eroanmam losık.aufsaker auytput,s? 

FE a.6b,.2,.ertyelltskoolean?: 

# 

kanmt nicht.wern E mıcht kommt 

und CE Koanmen Seneinzaäm oder san mıcht 
k:ammt ser Co Mmicht kommt und umgekehrt 


LmT 


E50 
2-7 nz 
for »3=true daunto Falze cd 


+cor Kbıetrue doupmto Falze da 


in 


I 
ıı 


+op cy=etrge dewntao Falze da 
I-I@ au 2) 
Fyelltssinat kilspot 2. and kzez) zarı) alba: 


r- 
fr er+tyuellt then urıtelnoa,.kb.ı. 
2) 


Pe 


oe ıdı 


she 


rc, 


a = 
Do 


a a 


, 
On 


=“ 


- 
| 


Co 


. 
Ce n 


De N) 


a- 
Ev 


v-, -—_ nen: 
ED an a9 Be BEE Ba un BT BE LT Er 


De) 


g 
5 


x) 


x 


Ze Dion) hin den nd hi ii MER EN Bis m ed di iu ie 
N 


m NIE GG ke 5 PD ii Ei in mn nn a 


m 
= 


FEN LOGIERUFGREE 

FRIHTCHREE 1) 

FEM A KOMMT HICHT.MEHNH B NICHT KOMMT 

FEM E UNE © KOMMEH GEMEINSAM ODER GAR MICHT 
FEM A KOMMT.MEHH © HICHT KOMMT UHO UMGEKEHRT 
For A=1 Tü -1 STEPR —- 
Für E=1 TO -1 STEPR - 
FÜR C=1 TO -1 STEF — 
F=-E:0=-R 

GOZUE Z7GsREM IMPLIERATIOH 
HI=WM:F=SE:0=SG 

GOSUB dZOsREM HEOUTYALEHNZ 
F=WH1 :0=) 
GOSUB Sa. 
YHz=W:F=Asdel 
GOZUE 4dEO:FKEM EAKLUSTIY-ÜODER 
P=H&:0=hl 

GOSUEB SagsrREM UHD 

IF YH=1 THEH PRIHTA,.E.C 

HE=T C 

HE=T E 

HE=T A 

EHL: 


fü fa fe 


M;FEEM UND 


„as E58 


REM IMPLIKATICHN 

IF Fi=0 THEN W=l:soTc 338 
M=-1 

RETURN 

REM AEBUIYALEHE 
WSSGNEPFO) 

RETURN 

REM ExKLUSIV-ODER 
W=-2GNEP#D) 

RETURH 

PEM UHE LOGIK —-AUFGSBEBE 
MESGHLP+D-1 3 

FETURM 


. 
= 

i 
„> 
pi 
[war 
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Berechnung von Determinanten 


von Karl Achilles 


] AUFGABENSTELLUNG 


Jeder quadratischen Matrix läßt sich eindeutig eine Zahl zu- 
ordnen, welche man als Determinante bezeichnet. Für eine 
(2,2)-Matrix läßt sich die Determinante auf einfache Weise 


berechnen. Die allgemeine Darstellung einer solchen Matrix ist 


a a 
= 11 12 1) 


Ay, Ay 


Für die Determinante |A| gilt dann: 


Al =a,, aa, 0 a," Ay, (2) 


Zur Berechnung einer n-reihigen Determinante (n >? 2) gibt es 
verschiedene Methoden, insbesondere die Entwicklung der Deter- 
minante nach Zeilen bzw. Spalten. Für die numerische Berech- 


nung mit Hilfe von Computern eignet sich jedoch eher ein Algo- 


2 BESCHREIBUNG DES LÖSUNGSWEGES 


Mittels Verdichtung wird eine Determinante n-ter Ordnung auf 
eine solche (n-1)-ter Ordnung reduziert. Durch sukzessive Re- 
duktion ergibt sich schließlich eine Determinante 2. Ordnung, 
welche sich nach (2) berechnen läßt. 

Die Methode der Verdichtung sei im folgenden erläutert. 

Eine n-reihige Determinante resultiert aus einer (n,n)-Matrix 
und hat folgende allgemeine Form: 


a a a ... a 


11 %ıa ®ıa ın 
aa, Ay2 Aa3 *** Ay 

A| = Ay] Ay2 9 rer A, (3) 
Ani a2 .. ..o aan 
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Voraussetzung ist, daß das Element a,;, von Null verschieden 
ist. Falls aber a,” O, so wird irgendein von Null verschie- 
denes Element ai durch Vertauschen von Spalte 1 mit Spalte i 
an die erste Stelle gebracht. Dadurch ändert sich das Vorzei- 
chen der Determinante, aber nicht ihr Betrag. Falls in irgend- 
einer Zeile (oder Spalte) alle Elemente gleich Null sind, so 
hat die Determinante den Wert Null! 


Ist die oben genannte Bedingung erfüllt, so werden folgende 
Operationen vorgenommen: 


- Die Zeilen 2 bis n werden mit a, multipliziert. Dadurch 
ergibt sich der an" !_Fache Wert der Determinante. 


- Zeile 2 -a - Zeile 1 


21 
Zeile 3 - az] - Zeile 1 
Zeile n - a1 « Zeile 1 


Man erhält folgendes: 


a a .oeeo. a 


11 12 in 
„n-1 a! _ OÖ A), ... Ay _; A), ... Ayn 
11 11 
0 Anz An Anz ._ An 
Daraus folgt schließlich: 
A), ... Ay 
_. _2-n . . 

al=a, |: : (4) 

Anz An 


Die Determinante, bestehend aus den A,yı ist nur noch von 
(n-1)-ter Ordnung. Für die A,x gilt: 


(i,k von 2 bis n) (5) 





Auf diese Weise läßt sich eine Determinante n-ter Ordnung auf 
eine solche der Ordnung 2 zurückführen. 
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3 VERBALE NOTATION DES ALGORITHMUS 


ANFANG 
DET:=1 
SOLANGE ZEILENZAHL > 2 TUE ANFANG 
FALLS A(1,1)=®, DANN ANFANG 
SUCHE GRÖSSTES ELEMENT DER 1. ZEILE (PIVOT) 
FALLS ALLE ELEMENTE DER 1..ZEILE NULL, DANN 
AUSGABE (DET = ®) ENDE 


VERTAUSCHE SPALTE DES PIVOTS MIT 1. SPALTE 
DET:= -DET 
ENDE 


FÜR ZEILE :=2 BIS ZEILENZAHL TUE ANFANG 
FÜR SPALTE :=2 BIS ZEILENZAHL TUE 
A(ZEILE,SPALTE) := A(1,1)*A(ZEILE,SPALTE) - 
A(ZEILE,1)*A(1,SPALTE) 
ENDE 
DET:= 
FÜR ZEILE:=2 BIS ZEILENZAHL TUE ANFANG 
FÜR SPALTE :=2 BIS ZEILENZAHL TUE 
A(ZEILE-1,SPALTE-1) := A(ZEILE,SPALTE) 


Der*a (1,1) (2"ZEILENZAHL) 


ENDE 
ZEILENZAHL:= ZEILENZAHL-1 
ENDE 
DET:= DET*(A(1,1)*A(2,2) - A(1,2)*A(2,1)) 
AUSGABE DET 
ENDE 
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4 PROGRAMMLISTINGS (BASIC und Pascal) 


BASIC-Programm Determinante 


ULIST 


1000 REM %* BERECHNUNG EINER N- 
REIHIGEN DETERMINANTE % 
1010 REM == 


1020 REM KARL ACHILLES 2.2.82 
1030 REM —- m nn 


1040 HOME =» INVERSE 

1050  FRINT "DETERMINANTE": NORMAL 
: VTAR (DD 

1060  INFUT "GEBEN SIE DIE ZEILEN 
ZAHL AN: ":N 

1070 IF N 2 THEN FRINT "N MUS 
5 GROESSER ALS 1 SEIN !": GOTO 
1060 

1080 DIM N,MA(N,N) „HILF (N,ND 

1090 FRINT =: FRINT "UND NUN DIE 
MATRIXELEMENTE (ZEILENWEISE) 
"= FRINT : FRINT 

1100 GOSUR 1240: REM EINGABE 

1110 HOME : VTARB (10): HTAB (10) 


1120 FRINT "EITTE WARTEN „.." 
1130 GOSUR 1340: REM BERECHNUNG 


1140 HOME 

1150 FRINT "DIE DETERMINANTE HAT 
DEN WERT "DET 

1160 END : REM ENDE HAUFTFROGRA 


MM 
1180 REM 


1190 REM %x%kxk UNTERFROGRAMME %%% 
1200 REM ------ nn 


1210 REM 

1220 REM 

1220 REM nenn 

1240 REM EINGABE DER MATRIXELEM 
ENTE 


1250 FOR IT = 1 TON 
1260 FOR J= 1 TON 


1270 PRINT "ACT". = "3 INFUT 
MACI,D 

1280 NEXT J,1l 

1290 RETURN : REM ENDE EINGABE 


1200. REM nn nn 


1310 REM 
1320 REM 


1370 REM mn nn 


1340 


1390 


1360 


1.370 


1380 


1390 
1400 


1410 
1420 


1430 
1440 


1490 


1440 
1470 


1480 
14970 
15300 
1510 
1520 
15.20 


1540 


1550 


1560 


1570 
1580 


1570 
1,00 
IA) 
16.0 
15.0 


1640 
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REM BERECHNUNG DER DETERMI 
NANTE 
HIDET = 1:ZZ = N 

IF ZZ = 2 THEN 1620: REM R 
ERECHNE Z2-REIHIGE DETERMINAN 
TE ELEMENTAR 

IF MA1,1) © 50 THEN 1520 


REM SUCHE FIVOTELEMENT (1. 
ZEILE) 


FÜR EL = 2 TO ZZ 

IF ABS (MACG,KID > ARS (M 
Atıl,E — 129) THEN FIVOT = MA( 
1,EI:MERE = HE 


NEXT E 

IF FIVOT = OÖ THEN DET = ©: 
RETURN 

REM 


FÜR LAUF = 1 TO ZZ: REM VE 
RTAUSCHE SFALTE 1 MIT SPALTE 
MERE. 
HILF (LAUF. MERK) = MACLAUF,ME 
RE.) 
MACLAUF,MERE) = MACLAUF, DD 
MAcLAUF, ID = HILF «LAUF, MERE) 


NEXT LAUF 


REM 

HIDET = - HIDET: REM ENDE 
FIVOTSUCHE 

REM 


FÜR ZE = 2 70 22 
FÜR SF = 2 TO ZZ 
MACZE,SF) = MAt1,13 % MACZE, 
SF) — MACZE.1? *% MAt1,Sr) 
NEXT SF.ZE 
HOCH = SGN (MA(1.,Ä1?) X ABS 
(MAG, = (ZZ - DD 
HIDET = HIDET 7 HOCH 
FÜR ZE = 2 TO ZZ: FOR SF = 
= TO ZZ 
MAIZE — 1,5SF - 12 = Mätlk,sh 
> 
NEXT SF,ZE 
22 — ZZ - 1: GOTO 14H 
DET - HIDET X MAL, 12D x Mt 
Zug = MAC1L.DI KK MAC, 1) 
RETURN = REM ENDE DETERMIN 
ANTENBERECHNLUNG 
REM un mn nn nn an nn 
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Beispiel 1 (BASIC) 


ÜRUN 
DETERMINANTE 
GEBEN SIE DIE ZEILENZAHL AN: > 


LIND NUN DIE MATRIXELEMENTE (ZEILENWEISE>» 
At1,12)= 


No 


A1,D)= 


D 
ul ua} 
ab 

n 
00 

! 


A222 
u 
Az, = 
v1 
AU3,1) 
AG 
Alı,2)= 
1 
A(},M= 


me 
rg 


EITTE WARTEN ... 
DIE DETERMINANTE HAT DEN WERT 7 


Beispiel 2 (BASIC) 


URLUN 
DETERMINANTE 
GEBEN SIE DIE ZEILENZAHL AN: 


UND NUN DIE MATRIXELEMENTE (ZETILENWEISE) 


All,l)= 
0 
All,D= 
20 
A0,D= 
0 
A2,l)= 
> 


A(2,2)= 


BITTE WARTEN ... 
DIE DETERMINANTE HAT DEN WERT 0 
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Pascal-Programm Determinante 


PRISRAR GE TI RMINANTE: 


ERNST MX 10% (& MAaxımMAlLE ZEN.EN- UND SFALIENZÄHL X) 
TYFE MATRIX = ARRaAYal. „MAX, 1..MAXu OF REAL: 


VAR I : INTEGEF: 
NATA s MATRIX: 


FRAPERURE 116261 OL, ME INTEGER: VAR MAalsMalhl da 
VER <EILE,SFALTE : INTEBER;: 


BEGIN 
FOR ZEILE: =1 ID L DO BEGSIN 
FÜR SPALTE: =1 TOM DO BEGIN 
WRITEU'MATC „ZEILE, ’.".SFALTE, "oı='oa 
REAL CMATAZEILE,.5SFAITEU) 
END; (iX OF SFALIE X 
END; (* OF ZEILE # 
WRITELN:WRITEL.N;G 
END; (x OF LTIESMAT 


FUNCTION DET (N: INTEGERZ VAR MAT=MATRIKI:REAL: 


VAR FIVOT,F.,O.HÜCH,HILFSDET : REäl: 
ZEILENZAHL. ZEILE, SFALIE,. MERKE, LÄUF.E,.L : INIEGER:; 
HILF : MATRIX: 


BEGIN 
HILFSDET: =1; 
ZEILENZAHL:=NG 


(X UMFORMUNG DER DETERMINANTE MITTELS "VERDICHIUNG' X 
WHILE ZEILENZAHL.-2 DO BEGIN 
IF MATAL,10=0 THEN BEGIN (x FIVOTSLICHE X 
FIVOT:=0; 
FOR K:=2 TO ZEILENZAHL DO 
IF ABSCMATAL.HEUI ABS (MATAI,K-10) THEN BEGIN 
FIVOT:=MATA1L.KO: 
MERKE: =E 
END; (x OF MAXELEMENT X) 
IF FIVOT=0 THEN BEGIN 
DET: =03 
EXIT(DEND: 
END; 
FOR LAUF:=1 TO ZEILENZAHL DU BEGIN 
‘%* VERTAUSCHE 5FALTE 1 MIT SPALTE MEREE % 
HILFAL ALF, MERH- Eues =MAlALALUF, MERKEus 
MATaLAUF, MERKEEÜ: =MATALAUF, 10: 
MATAL AUF, lu: =HILFALAUF, MERE Eu 


END; (X UF VERTAUSCHE 3% 
HILFSDET: =-HILFSDET; 
END; (x OF FIVOTSLCHE X 


63 


Determinanten 





(X VERDICHTUNG 

F:=MATA1,10; 

FOR ZEILE:=2 TO ZEILENZÄHL DO BEGIN 
Ü:=MATAZEILE, 10: 
FOR SFALTE:=2 TOD ZEILENZAHL DO 


MATAZEILE,.SFALTEÜ: =FXMATAZEILE, SFAL.TEU-OKMATAL,SFALTEU 


END; (X OF ZEILE 

(%X ENDE VERDICHTUNG 
HOCH: =MATAL, 103 
L:=ZEILENZAHL-2%5 


WHILE L>1 DO BEGIN (% BERECHNE MAT HOCH (ZEILENZAÄHL-2I 8% 


HOCH: =HOCHXMATAIL. 1Us 
L:6-1 


END; (* OR MAT HOCH 
HILFSDET: =-HILFSDBET/HOCH; 
FÜR ZEILE:=Z TO ZEILENZAHL DO BEGIN 
(% REDUZIERUNG DER MATRIXDIMENSION % 
FOR SFALTE:=2 TO ZEILENZAHL DO 
MATAZEILE-1,S5FALTE-1U:=MATAZEILE,.SFALTEU 
END; (X DF REDUZIERUNG % 
ZEILENZAHL. : =ZETLENZAHL -1 
END: (* OF (ZEILENZAHL>2I OR (NOT AUSGANG) 
DET: =HILFSDETK (MATAL,. LUXMATAZ, ZU-MATA1, ZUXMATAZ, 109; 
END; (x OF DET X 


BEGTN (X HAUFTFROGRAMM %) 


WRITE (CHR (12); (% LOESCHT BILDSCHIRM 
WRITELN(C"DETERMINANTE"3G 

WRITELN(’"==s=s=er2m 2222: }) 4 

WRITELN:WRITELN: 

WRITE’ ZEILENZAHL 7'235 

FEADCIIDG: 

L1I1ESMAT(IL, I1,MATAI: 

WRITELNC'"DETERMINÄNTE — °„DEI(IL,MATAI:S 12:6) 


END. 


5 LITERATUR 


x) 


Aitken, A.C.: Determinanten und Matrizen. Mannheim: 


Bibliographisches Institut 1969 
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Polynomerzeuger 


Polynomerzeuger 


von Karl Achilles 


1 AUFGABENSTELLUNG 


Ein Polynom P n-ten Grades habe n reelle Nullstellen, welche 
bekannt seien. Gesucht sind die Koeffizienten a,(i = 1 bis n) 


des normierten Polynons. 


2 BESCHREIBUNG DES LÖSUNGSWEGES 


2.1 Beispiel 

-2; 1; 3 seien die Nullstellen. 

Das normierte Polynom hat dann die Linearfaktorzerlegung 
P = (x+2) (x-1) (x-3) = x°-2x7-5x+6. 


Ergebnis a,=1 a,=2 a,=-5 a=6 


2.2 Verallgemeinerung 


Xi Kgi er 5 X_45 %. seien die n Nullstellen eines Polynoms 
n-ten Grades. 
Es gilt dann: 


x + ... ta x2ta xta_. 
- 2 1 (6) 


Nach dem Wurzelsatz von VIETA sind die a, mit den x, folgen- 


dermaßen verknüpft: 


a = 1 
n 
= - (x, tx_+ ... +x 
An-1 172 n) 
a = x X_+x\X-+X_X-+ ... +X 


n-2 12 1,3 72% 3 n-1“n 


= (RK RR 3X RX gr ... +X 


[e)) 
) 


n-2*n-1°n) 


n 
o (-1) "XıKgXZ" oo. X 4%, 


n 
I 


© Friedr. Vieweg & Sohn Verlagsgesellschaft mbH, Braunschweig 1983 65 


Polynomerzeuger 





3  PROGRAMMBESCHREIBUNG 


wie das folgende Struktogramm verdeutlicht, kann man für den 


VIETAschen Wurzelsatz einen einfachen Algorithmus angeben. 


3.1 Struktogramm 


Eingabe des Polynomgrades n und der Nullstellen X, 






(für i=1 bis n) 










Ausgabe der Koeffizienten a, (für i=@ bis n) 


Daß dieser Algorithmus das Problem löst, wird im folgenden 





exemplarisch gezeigt. In der Tabelle wird der Berechnungsteil 


aus dem Struktogramm (für die a, für i=® bis n-1) durchge- 


spielt. 

zugrunde liegt das Beispiel: n = 3 
x = 
X, = 2 
X, = 3 


Zu Beginn der Rechnung gilt: a,= und a,=® (i=® bis 2) 
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rn rar ar 


Tabelle (Belegungsplan): 





a7 
ao 

3.2 Programmlistings 

urR#U 

Wl_IST 

IiH00 REM FÜLYNOMERZEUGER 

1010 REM KKEKKKKRKKEKKRKKKKK 

1020 REM KARL ACHILLES 14.1.8232 

xx% 

1030 REM 

1040 HOME 

1650 INVERSE : FRINT "FOLYNOMERZ 
EUGER": NORMAL 

i0O&0 VTAR (5) 

1070 REM 

1080 REM EINGABE 

1100 INFUT "GEBEN SIE DEN GRAD D 
ES FOLYNOMS AN „. "ıN 

1105  DIM X{N),.AN 

1110 FRINT : FRINT "GEBEN SIE NU 
N DIE NULLSTELLEN AN :" 

1120 FOR IT = 1 TON 

1120  FRINT I".NULLSTELLE .„. "; 

1132 INFUT X(D 

13135 NEXT I 

1127 REM ENDE EINGABE 


1. 


Nullstelle 


wird verwendet 


2. Nullstelle 
wird verwendet 


3. Nullstelle 


wird verwendet 


11.:8 
1139 


1140 
1150 
1155 
1157 
1150 
1170 
1190 


1190 
1195 
1177 
1200 


1205 


1210 
12.0 
12.30 


1.50 


REM 

REM BERECHNUNG DER KOEFFI 
ZIENTEN 

LETAM = 1 

FÜR T=0TON -— 1 

LET AUT) = © 

NEXT I 

FÜR E = N TO 1 SIEFR - 1 
FOR IT = N TO 1 S5STEF — 1 
LET A{N- DD = AN - Tı ca 
A(iN- I + 1 X X{N- K + 19) 
NEXT 1,E 

REM ENDE BEREIHNLUNG 

REM 

REM AUSGABE DER AD 

FRINT =: FRINT = FRINT "FNEF 
FTIZIENTEN DES FÜLYNOMS": FRINI 
FAR I =N TO © STEF —- 1 
PRINT "At"1l") "Aıl) 

NEXT ] 

END 
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: Pokal FO VNOMERZEIGERS 
(al MAX IMALGRAD 0 (X WILLEUERLICH GEWAEHLT %) 


12 SReaß, Ta HS INTEGER;: 
HILF ; REAL: 
NULLSTELLE 3 HRRAYAl..MAXIMALGRADU UF REAL; 
FORFFIZTENT :5RRAYAO. „MAXIMALGRABU OF REAL; 
SEUON 
(X UFBERSCHRIFI X) 
RTTE CHF EIEIDG (% SAEUBERT BILDSCHIRM 
STTELN «FOL YNOMERZEUGER’); 
URTTELN "een een) 
WRITELN:WRITELN: 


(% EINLABE % 

WRITE C"RTTTE FOLYNOMGRAD EINGEBEN: "34 
F}ADILN (GFADIE WRITELN: 
We TTELH CHLTTE NUNLSTELLEN EINGEBEN’; 
KHLTELN 6° mem nn nn — wu 
KIELTELN: 
FIR HE =» 1 TO GRAD DO 
WLOEN 

WRITE C’NULLSTELLE NR.’,K,’ = °9; 

READLN (NULLSTELLEAKOIG 
FiID: 


(X FOEFFIZIENTENANFANGSWERTE X 
Für T e= 6 TO GRAD-1 DO 
FOEFFIZIENTAlU x: 05 
IEFFIZITENTAGRADUL := 1; 


ı%X KOEFFIZIENTENBERECHNUNG X 
FÜR Fo+= GRAD DOWITO 1 DO 
HHGIN 
FÜR I := GRAD DIIWNTO 1 DO 
REGIN 
HILF := NULLSTELLEAGRAD-K+lüs 
FOEFFIZIENTÄGRAD--Iü 2:= KOEFFIZIENTABGRAD- Iu-KOEFFIZIENTAGRAD- 1 +IUXHILF 
END; 
END: (% ENDE KOEFFIAZIENTENBERECHNLUNG % 


ı* KDOEFFIZIENTEN AUSDRUCKEN % 
FOUR I :s= 1 TO 3 DO WRITELN: 
FOR I ;= GRAD DOWHTD 0 DO 


WRITELN C’KOEFFIZIENT A(’,1,') "„KOEFFIZIENTAIU)G 


ı ND. 
4 WEITERE BEISPIELE 


4.1 Nullstellen -4 -2 -1 1 2 


Koeffizienten a,=1 a,=4 a,=-5 a,=-29 


4.2 Nullstellen 0.5 -1.5 -6.1 0.01 ) 
Koeffizienten a,=1 a,=7.99 a,=5.279 


s8 a,=-4.6285 a,=9.4575 a,=® 


Größter gemeinsamer Teiler 


Größter gemeinsamer Teiler (ggT) 


von Karl Achilles 


1  AUFGABENSTELLUNG 


Von zwei natürlichen Zahlen a,b soll der ggT bestimmt werden 


nach dem Euklidischen Divisionsalgorithmus. 


2 LÖSUNGSWEG 


Der Euklidische Divisionsalgorithmus basiert auf einer fort- 
laufenden Teilung mit Rest. Er terminiert, wenn der Rest Null 


geworden ist. Dies sei an einem Beispiel veranschaulicht: 
Zu bestimmen sei ggT(437;69): 

437 : 69 = 6 Rest 23 

69 : 


IN 
II 
Il 


3 Rest_Q_ 


Also ist ggT(437;69) = 23 


3  PROGRAMMBESCHREIBUNG 


3.1 Struktogramm 


Lies (A,B) 


Wiederhole 












bis R=® 


Schreibe (A) 
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Größter gemeinsamer Teiler 


I 
3.2 Programmlistings 


Größter gemeinsamer Teiler: BASIC-Programm 


ULIST 

100 REM x 66T % 
110 REM 

120 HOME 


130 INVERSE »: FRINT "66T": NORMAL 


140 VTAR (59 
150 PRINT "GEBEN SIE BITTE ZWEI 
NATUERLICHE ZAHLEN A,RB EIN : 


1860 INPUT "A "aA 

170  INFUT "RB ": 

180 LET GSET = A: LETV=R 

1970 REM 

200 REM BERECHNE G6GT 

210 LET REST = 66T - v% INT (6 
GT /W 

220 LET G6T = V 

2350 LET V = REST 


|| 


240 1IF REST > OÖ THEN 210 

250 REM 

260 REM AUSGABE 

270 HOME 

280 FRINT "G6T("A"S"R"Y) = "66T 
>00 END 


Beispiele in BASIC: 


URUN 

GGT 

GEBEN SIE RITTE ZWEI NATUERLICHE ZAHLEN A,RB EIN 
Ah = 75 

B = 15 

66T (75:15) = 15 


URUN 

GGT 

GEBEN SIE BITTE ZWEI NATUERLICHE ZAHLEN A,B EIN :; 
A = 87654321 

B = 12345678 

66T (876543215 12345678) = 9 
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Größter gemeinsamer Teiler: Pascal-Programm 


FROGRAM GGTEILER; 


VAR ZAHLI. ZAHLZ, REST, G6GT, V :INTEGER; 


BEGIN (% HAUFTFROGRAMM 


(X EINGABE % 

WRITELN C’BITTE 2 NATUERLICHE ZAHLEN EINGEBEN’)D3 
READ (ZAHL.1,ZAHLED ; 

GGT := ZAHL.1ls V = ZÄHLE; 

(k BERECHNUNG DES 66T 


REFEAT 
REST := GGT MOD V; 
GGT ı= Vi V s= REST 
UNTIL REST = 0: 


(% AUSDRÜCKEN DES G6T % 

WRITELN; 

WRITE €°’GGT (°,ZAHL1,’3"„ZAHL2,’) = °,667; 
END. 


3.3 Bemerkungen zur Programmierung 


Für die Division mit ganzzahligem Rest gibt es in Pascal die 
Modulo-Funktion (REST = GGT MOD V). Zum Beispiel ergibt sich 
für GGT=17 und V=5 der ganzzahlige Divisionsrest 2; also ist 
17 MOD 5 = 2. 

Die Modulo-Funktion muß in BASIC erst simuliert werden. Unter 
Verwendung der BASIC-Funktion INT (Ganzzahliger Teil einer 
zahl) kann man folgende Beziehung für den ganzzahligen Divi- 


sionsrest aufstellen: 
REST = GGT - V-INT(GGT/V) (3.31) 


Für GGT=17 und V=5 ist dann INT(GGT/V)=3 und REST=17-5-3=2. 
Die Beziehung (3.31) ersetzt also die Formel REST=GGT MOD V. 
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Monte-Carlo-Pi 


von Karl Achilles 





1 AUFGABENSTELLUNG 


Die Kreiszahl Pi soll mit Hilfe eines Zufallsexperiments ge- 


schätzt werden (Monte-Carlo-Methode). 


2  LÖSUNGSWEG 


Dazu läßt man N Regentropfen auf 

ein Quadrat der Seitenlänge 1 

fallen (siehe Skizze) und zählt 

die Anzahl Z der in das Innere | 
des Viertelkreises fallenden 

Tropfen. 

Der Flächeninhalt des Viertel- 

kreises mit dem Radius | 

entspricht Pi/4. 

Der Bruch Z/N stellt den Anteil 

der in den Viertelkreis fallenden Tropfen dar und ist daher 


eine Näherung für Pi/4. Es ergibt sich folgende Beziehung: 


. 4-2 
Pi» —m 


Der Simulation des Zufallsexperiments liegt folgende Idee zu- 
grunde: Man erzeugt mit Hilfe eines Zufallsgenerators hinter- 


einander zwei Zufallszahlen x,y im Bereich zwischen ® und 1. 


Denkt man sich den Ursprung eines Koordinatensystems in die 
linke untere Ecke des Quadrats gelegt, so wird durch die Ko- 
ordinaten x,y ein Punkt innerhalb des Quadrats beschrieben. 
Dieser Punkt liegt genau dann innerhalb des Viertelkreises, 
wenn die Bedingung x2 + y2 < 1 erfüllt ist. Alle Punkte mit 
den Koordinaten x,y, die auf dem Kreis liegen, erfüllen be- 
kanntlich die Kreisgleichung x? + y2 = 1. 

Es sollte nicht unerwähnt bleiben, daß die Monte-Carlo-Metho- 
de sehr schlecht konvergiert und daß unbedingt ein guter Zu- 


fallsgenerator verwendet werden sollte! 
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3  PROGRAMMBESCHREIBUNG 


3.1 Struktogramm 


Für LAUF von 1 bis N tue 











Erzeuge Zufallszahl x 
Erzeuge Zufallszahl y 


m 7 


Dnmmen nem | oo 
PI &— 4 -ZÄHLER/N 


Schreibe (PI) 


3.2 Programmlistings 
Monte-Carlo-Pi: BASIC-Programm 








Monte-Carlo-Pi 


UPR#O 
LIST 
100 REM % MONTE-CARLO-FI X 
110 REM 
120 HOME 
120 INVERSE : FRINT "MONTE-CARLO 
-FI": NORMAL . . . 
140 PRINT z PRINT Beispiele in BASIC: 
150 INFUT "WIEVIELE VERSUCHE ";N 
URUN 
1&0 LET ZAEHLER = 0 MONTE-CARLO-FI 
170 FOR LAUF = 1 TON 
180 LETX = RND (i): LETY= RND 
(1) WIEVIELE VERSUCHE 10 
190 IFXXX-+-YXxY = 1 THEN ZA NAEHERUNG FUER FI = =.2 
EHLER = ZAEHLER + i 
200 NEXT LAUF 
210 LET FI = 4 *% ZAEHLER / N 
220 FRINT "NAEHERUNG FUER FI = " 
rl ÜRLIN 


AO 


END 


MONTE-CARLO-FI 


WIEVIELE VERSUCHE 1000 
NAEHERUNG FUER FI = 3.128 
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Monte-Carlo-Pi: Pascal-Programm 


FROGRAM MONTEPI ; 


USES AFFLESTUFF; (%k ENTHALT ZUFALLSGENERATOR %) 
VAR ZAEHLER,LAUF,MAX : INTEGER3 

X.Y,PI : REAL3 
BEGIN 


WRITE(C’"WIEVIELE REGENTRÜFFEN 7 "93; 
READLN (MAX; 
ZAEHLER: =0; 
RANDOMIZE; 
FOR LAUF:=1 TOD MAX 
DO REGIN 
(k ERZEUGE ZUFALLSZÄHLEN % 
X: =RANDOM/FE27673; 
: =KANDOM/S32767; 
IF X%XX+YxY<1 THEN ZAEHLER: =ZAEHLER+1; 


END; 
Fi:=4XZAEHLER/MAX; 
WRITELN;: 
WRITELNÜ’NAHERUNG FÜR FI = °,FD; 
WRITELN; 
WRITELN(MAX„" ZUFÄLLSVERSUÜCHE”); 


END. 


3.3 Bemerkungen zur Programmierung 


Der wesentliche Unterschied in den Programmen (BASIC, Pascal) 
besteht in der Erzeugung der Zufallszahlen: 

In Pascal wird durch die Funktion RANDOM eine ganzzahlige Zu- 
fallszahl zwischen ® und 32767 erzeugt, welche mittels Divi- 

sion durch 32767 auf das Intervall [9;1] transformiert werden 
muß. In BASIC steht durch die Funktion RND(1) sofort der ge- 


wünschte Bereich zur Verfügung. 


4 LITERATUR 


Enge |], A.: Elementarmathematik vom algorithmischen Stand- 
punkt, Klett 1977 
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D’Hondtsches Höchstzahlverfahren 


D’Hondtsches Höchstzahlverfahren 
von Karl Achilles 


1] AUFGABENSTELLUNG 


Nach einer Wahl sollen M Mandate auf N Parteien nach dem 
D'Hondtschen Verfahren verteilt werden. Dabei werden die 
Zweitstimmenzahlen jeder Partei nacheinander durch 1,2,3,... 
dividiert und anschließend die Mandate in der Reihenfolge der 
Größe der anfallenden Quotienten vergeben. Hat eine Partei we- 
niger als 5% aller Stimmen bekommen, so erhält sie kein Man- 
dat. 


2 LÖSUNGSWEG 


Für die Zweitstimmenzahlen der N Parteien wird ein Feld Zur 
Zur 231 ... Zur bereitgestellt. Entsprechendes gilt für die 
Quotienten Q,; und die Divisoren (Teiler) T, (i=1 bis N). 
Zunächst werden die Zweitstimmenzahlen den Quotienten zugeord- 
net, wobei alle Teiler den Wert 1 erhalten. Es gilt also: 

Q; = zZ; für i von 1 bis N. Befindet sich das Maximum aller 
Quotienten an der k-ten Stelle, so wird der k-ten Partei das 
erste Mandat zugeteilt. Der Quotient Q, wird gestrichen und 
durch einen neuen Quotienten ersetzt, welcher sich aus 
Qu=2,/T, mit T,=2 ergibt. Somit hat sich für die k-te Partei 
der Teiler um 1 erhöht, nachdem ihr ein Mandat zugeteilt wor- 
den war. 

Nun wird das nächste Maximum der Q, (i=1 bis N) gesucht. Nach- 
dem ein Mandat zugeteilt wurde, muß immer der entsprechende 
Teiler um 1 erhöht und der neue Quotient gebildet werden. 
Dieser ersetzt dann den vorherigen Quotienten. 

Nach jeder Zuteilung wird der Mandatszähler (der alle Mandate 
zählt) erhöht. Das Programm bricht mit der Berechnung ab, 

wenn gilt: Mandatszähler = M. 
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D’Hondtsches Höchstzahlverfahren 


3  PROGRAMMBESCHREIBUNG 


3.1 Struktogramm 


Anfang 


Lies (PARTEIENANZAHL, MANDATSANZAHL) 


Für I von 1 bis PARTEIENANZAHL tue 








Lies (Parteiname (I)) 
Lies (ZWEITSTIMMENZAHL(I)) 








Berechne die GESAMTSTIMMENZAHL 
Berücksichtige die 5%-Klausel 






Für I von 1 bis PARTEIENANZAHL tue 











MANDATZAHL (I) 
TEILER(I) <— 1 
QUOTIENT (I) <— ZWEITSTIMMENZAHL (I) 


<— 9 













MANDATZUTEILUNG <— ® 


MANDATZUTEILUNG <— MANDATZUTEILUNG+]I 
MAXIMUM <— OQUOTIENT (1) 
Für I von 2 bis PARTEIENANZAHL tue 


wiederhole 










QUOTIENT (I) > MAXIMUM 


Anxımm <— ovorzenmn | 7 


Erhöhe MANDATZAHL der Partei mit MAXIMUM 
Erhöhe TEILER " " " " 
Berechne QUOTIENT " " " " 











bis MANDATZUTEILUNG 


Für I von 1 bis PARTEIENANZAHL tue 







Schreibe (PARTEINAME(I)) 
Schreibe (MANDATZAHL(I)) 
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3.2 Programmlistings 
D'Hondtsches Verfahren: BASIC-Programm 


üc 

ULIST 

1000 REM KK -—- D’HONDT -— xık 

1010 REM 

1020 REM KARL ACHILLES 9.11.81 

1050 REM 

1040 REM ---------- -- - ---- - 

1050 REM ZI) ... ZWEITSTIMMENZ 
AHLEN 

1040 REM N$(I)... PARTEINAMEN 

1070 REM O@(D ... BUOTIENTEN DE 
RI-TEN FARTEI 

1080 REM K(D .„.. TEILER FUER D 
IE I-TE FARTEI 

1090 REM M(T .„... ANZAHL DER MA 
NDATE FUER DIE I-TE FARTEI 

1100 REM ---------- - ------ 

1110 REM 

1120 REM ------------ - - --- 

11350 REM HAUFTFROGRAMM 

1140 REM --------- - - - -- 

1150 REM 

1140 HOME 

1170  FRINT "D’HONDTSCHES VERFAHR 
EN" 

118BÜ FRINT "==s==s=-=<-=- === mes =meoen 

1170 REM 

200  GOSUR 1290: REM EINGABE 

1210 GOSUR 1460: REM GESAMTSTIM 
MEN 

1220  GOSUR 1570: REM S5%-KLAUSEL 

1250  GOSUR 16560: REM INITIALISI 
ERUNG 

1240 GÜSUR 1740: REM BERECHNE M 
ANDATSVERTEILUNG 

1250 GOSUR 1940: REM AUSGABE 

1260 END 

270 REM ENDE HAUFTFROGRAMM 

1280 REM 

1290 REM 

1300 REM --- EINGABE --- 

1310 REM ---- nn 

13220 VTAR (5 

1220 INFUT "GEBEN SIE DIE ANZAHL 


1340 
1350 


1360 
1370 


1380 


1390 


1400 
1410 


14.30 
1440 


1450 
1460 
1470 
1490 
1490 


1500 
1510 


8) 


1 


{N 
rJ 


fr 


Do 


ee 

£n ca ch 

a u | 
x 


10) 


1560 
1570 


1580 
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DER FARTEIEN AN "N 

INFUT "GEBEN SIE DIE ANZAHL 

ALLER MANDATE AN":MA 

DIM Z(ND .„N$(ND .DCN).KAN) Mt 
N) 

VTAR (15) 

FRINT "LESEN SIE DIE NAMEN 
DER FARTEIEN EIN -" 

FRINT "ZUERST DIE FARTEI MI 
T DER HOECHSTEN STIMMENZAHL"” 
: FRINT 

FOR I = 1 TO Ns INFUT N$(D 
: NEXT I: FRINT 

HOME 

FRINT "GEBEN SIE DIE ZWEITS 
TIMMENZAHLEN DER FARTEIEN EI 
N :"s FRINT 

FOR I = 1 TO N: FRINT "ZWEI 
TSTIMMENZAHL DER FARTEI "N$( 
ID: INFUT Z(I)s NEXT I 

RETURN 
REM --- ENDE EINGABE --- 


REM 

REM --- BESAMTSTIMMEN --- 
REM ------- 4. 
LETG = 06 

FOR I = 
Z(Id: NEX 
HOME 
FRINT "DIE SUMME DER ZWEITS 
TIMMEN IST: "G 

FRINT "--------- 
RETURN 

REM 

REM --- ENDE GESAMTSTIMMEN 


_— 


| m 


REM 

REM --- BERUECKSICHTIGUNG 
DER S%A-ELAUSEL --- 

REM ------- nn 
LET LIMIT = .05 x6G 

FOR J= 1 TON: LET Z(W) = 
ZU X SGN ( INT (Z(I) / LI 
MIT): NEXT J 

RETURN 
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D'Hondtsches Verfahren: BASIC-Beispiel 


1620 REM 1860 LET SM = SM + 1: REM SM IS 

1650 REM --- ENDE S%-KLAUSEL == T ZAEHLER FUER DIE ZUGETEILT 
” EN MANDATE 

1640 REM 1870 IF SM = (MA) THEN GOTO 192 

1650 REM --- INITIALISIERUNG -- o 


_ 1880 REM ELSE ... ERHOEHE TEILE 
1660 REM een R,EILDE NEUEN QUOTIENTEN UND 
- DURCHLAUFE WEITER DIE SCHLE 


1670 LET SM = OÖ IFE 
1680 FOR I = 1 TON: LETKID = 1896 LET K(MERKE) = K(MERKE) + 1 
1: NEXT I 
1696 FOR I = 1 TON: LET BID = 1700 LET R(MERKE) = Z(MERKE) / K 
1700 "RETURN (NERKE) 
on 1910 G0TO 1770 
1710 REM --- ENDE INITIALISIERU 1920 RETURN 
1720 REM 1930 REM --- ENDE BERECHNUNG -- 
1720 REM 1940 REM 
1740 REM --- BERECHNE MANDATSVE 190 REM 
a > 1960 REM --- AUSGABE --- 
IT 1970 REM ---------- 
nr 1980 VTAR (10) 
1760 REM 1990 PRINT "DIE MANDATE VERTEILE 
1770 LET MERKE = 1 . 
1780 LET HILF = @(1) N nur 
- = tl 2000 PRINT "EINZELNEN FARTEIEN W 
1790 FOR I = 1 TON IE FOLGT :" 
1800 IF @(D < = HILF THEN G0DTO 010 PRINT "=eemenaanaenaneennnn 
1840: REM VERGLEICH DER @(I ET neeseseoeeeenene u 
zamzzmmzemmsmensoene ı FRINT 
) MIT @(1) © PRINT 
1810 REM ELSE ... VERTAUSCHE 2020 FOR I = 
la20 LET H Soc 2020 FOR I = 1 TON 
1920 LET MERKE = I 2030 PRINT TAB( 1DN$(D; TAB( 1 
2 = &)M(I); TAB( 21) "MANDATE" 
1840 NEXT I 2040 NEXT I 
1850 LET M(MERKE) = M(MERKE) + 1 5050 RETURN 
: REM TEILE MANDAT DER FART Sol ___ ___ 
EI NR.MERKE ZU 2060 REM ENDE AUSGABE 
ÜRUN 
D’HONDTSCHES VERFAHREN 
GEBEN SIE DIE ANZAHL DER FARTEIEN AN 3 
GEBEN SIE DIE ANZAHL ALLER MANDATE AN7B 
LESEN SIE DIE NAMEN DER FARTEIEN EIN - 
ZUERST DIE FARTEI MIT DER HÜECHSTEN STIMMENZAHL 
PEDU 
?SFD 
>FDF 
GEBEN SIE DIE ZWEITSTIMMENZAHLEN DER FARTEIEN EIN : 
ZWEITSTIMMENZAHL DER FARTEI CDU 
7987654 
ZWEITSTIMMENZAHL. DER FARTEI SFD 
7894566 
ZWEITSTIMMENZAHL DER PARTEI FDF 
?190106 
DIE SUMME DER ZWEITSTIMMEN IST: 2072326 
DIE MANDATE VERTEILEN SICH AUF DIE ED E NE 
N FARTEIEN WIE FOLGT : = = 
EINZELNE TEIEN 8 FDF 7 MANDATE 


Ionean Emma tem mann tetee Fitem snarı snae CpERe Amane HERE KENN sueen FröRS Meet Aähm Seste Sahne Sana MEER wiER mass Anane enane onane Mnue> His} anacs ment MiREE dasce vorm> Surme mann sadae Fidan msare munet Söre 
a een Suse Casmm Hrnan M00iE Senn ammen oanen mens Aanze Codes MIBER maNME aenen Onane Aare pitum FOHte aitEn name SAE Adacn Cases männl SEtie wesen enide 
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D'Hondtsches Verfahren: Pascal-Programm 


FRÜOGRAM DHONDT; 


(% DHONDTSCHES HOECHSTZAHLVERFÄAHREN %) 


CONST MAX = 105 (X MAXIMALE ANZAHL DER FARTEIEN % 


TYPE GANZBEREICH = ARRAYAl..MAXU OF INTEGER; 
REELLBEREICH = ARRAYAl..MAXU OF REAL; 
NAMENBEREICH = ARRAYAl..MAXU OF STRING: 
FARTEIEN = RECORD 

FNAME : NAMENBEREICH; 

STIMMENZAHL, 

TEILER, 

QUOTIENT : REELLBEREICH; 

MANDATZAHL : GANZBEREICH; 
END; (X OF FARTEIEN %) 

VAR FANZAHL, MANDATE, ZEIGER,LAUF : INTEGER; 
ALLE,LIMIT : REAL; 

FARTEI : PARTEIEN: 


FROCEDURE LOESCH; 
VAR HE : INTEGER; 
BEGIN 
FOR K:=1 TO 1000 DO RBEGIN END; 
WRITE (CHR (12993; 
END: (x OF LOESCHEN %) 


FROCEDURE EINGABE; 
BEGIN 
LOESCH; 
GOTOXY (10,109; 
WRITE (°’D’ ’HONDTSCHES VERFAHREN’); 
GOTOXY (10,12; 


LOESCH; 
GOTOXY(1,10); 
WRITEC"WIEVIELE FARTEIEN ?’); 
GOTOXY (30,10); 
READLN (FANZAHL> 5 
GOTOXY(1,1D3; 
WRITE("WIEVIELE MANDATE INSGESAMT 7? ”)5; 
GOTOXY(30,193; 
READLN (MANDATE) 3: 
END: (X OF EINGABE %) 


FROCEDURE LIES: 
FPROCEDURE LIESSTIMMEN; 
VAR I2 : INTEGER; 


BEGIN (%* LIESSTIMMEN %) 
WITH FARTEI 
DO BEGIN 
FOR I2:=1 TO PANZAHL DO 
BEGIN 
GOTOXY (0,12; 


WRITE (’ZWEITSTIMMENZAHL DER FARTEI ",„FNAMEAIZU,’ 


GOTOXY(35, 129; 
READ (STIMMENZAHLAI2UD: 
END: 
END; 
END; (* OF LIESSTIMMEN %) 


”) 


“ 
a 
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BEGIN (X LIES % 
LOESCH; 
FOR LAUF:=1 TO FANZAHL 
DO BEGIN 
GOTOXY(1,LAUFIS: 
WRITEC’NAME DER ",LAUF,’ „ FARTEI ? 93; 
GOTOXY(SO,LAUPDS 
READLN (FARTEI.FNAMEALAUFÜI: 
END; 
LOESCH; 
LIESSTIMMEN; 
END; (x OF LIES %) 


FROCEDURE ALLESTIMMEN; 
VAR IS : INTEGER; 
BEGIN 
LOESCH: 
ALLE: =0;3 
FOR I3:=1 TO FANZAHL 
DO ALLEs=ALLE+FARTEI.STIMMENZAHLAISU: 


WRITELN; 
WRITELN (’GESAMTSTIMMENZAHL = °,ALLEIS 
END; (X OF ALLESTIMMEN % 
FROCEDURE FUENFFROZENT; 
VAR I4 : INTEGER; 
BEGIN 


LIMIT:=0.0SXALLE; 
FOR I4:=1 TO FANZAHL 
DO IF FARTEI.STIMMENZAHLAI4U © LIMIT THEN FARTEI.STIMMENZAHLAI4U: =; 


END; (%X OF FUENFFRDOZENT % 
FROCEDURE INITIALISIERE: 

VAR IS : INTEGER; 

EBEGIN 


FOR 15:=1 TO FANZAHL 
DO WITH FARTEI 


DO BEGIN 
MANDATZAHLAISU:=0; 
TEILERAISU: =1; 
QUOTIENTAISU:=STIMMENZAHLAISU: 


END; 
END; (x OF INITIALISIERE % 
FROCEDURE BERECHNE ; 
VAR Iö, MANDATZUTEILUNG, ZEIGER : INTEGER; 
MAX IMUM : REAL; 
BEGIN 
WITH FARTEI 
DO BEGIN 
MANDATZUTEILUNG: =0;3 
REFEAT 


MANDATZUTEILUNG: =MANDATZUTEILUNG+135 
ZEIGER: =1; 
MAX IMUM: =QUOTIENTAIU; 
FOR Iö:=2 TO FANZAHL DO 
IF QUOTIENTAIGU > MAXIMUM THEN 
BEGIN 
ZEIGER:=I6; 
MAX IMUM: =QUDTIENTAISU 
END; (x OF IF WM 
MANDATZAHLAZEIGERÜ: =MANDATZAHLAZEIGERU+1; 
(% NEUES MANDAT ZUGETEILT % 
TEILERAZEIGERÜU: =TEILERAZEIGERÜU+135 
(% TEILER ERHOEHT %) 
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QUOTIENTAZEIGERU:=STIMMENZAHLAZEIGERUÜ/TEILERÄAZEIGERÜ; 
(%X NEUEN QUOTIENTEN BERECHNET % 
UNTIL MANDATZUTEILUNG = MANDATE; 


END; 

END; (k OF BERECHNE x) 
FROCEDURE AUSGABE; 

VAR 17 x: INTEGER; 

BEGIN 


WRITE <CHR (129); 
GOTOXY (5,0935 
WRITELN<(’MANDATSVERTEILUNG :°”2; 
GOTOXY (5.793 
WRITELN (’=s=====2=2=2=2m=2=2=2=m2m=mmm=? )3 
FOR 17:=1 TO FANZAHL 
DO BEGIN 
LAUF:=1I7+10; 
GOTOXY(1,LAUFPD; 
WRITE (FARTEI.FNAMEAI?UI:G 
GOTOXY(ZO,LAUF;S 
WRITE <(FARTEI. MANDATZAHLAIT7UIG 
GOTDXY (25,LAUFDG; 
WRITE (’MANDATE’I;: 
END; 
END; (kX OF AUSGABE 9% 


BEGIN (%k HAUFTFROGRAMM %) 
EINGABE; 
LIES; 
ALLESTIMMEN; 
LOESCH; 
GOTOXY (10,10); 
WRITEC"BITTE WARTEN „..’); 
FUENFFROZENT; 
INITIALISIERE; 
BERECHNE ; 
AUSGABE; 

END. 


3.3 Bemerkungen zur Programmierung 


Einer Partei, die weniger als 5% der Gesamtstimmenzahl er- 
reicht, wird die Zweitstimmenzahl ® zugewiesen. Dadurch erge- 
ben sich alle Quotienten zu Null, und die Partei erhält kein 
Mandat. 

In der Pascal-Version werden alle Daten einer Partei mit Hilfe 
des Datentyps RECORD (Verbund) zusammengefaßt. Dieser Typ er- 
möglicht den Zugriff auf sämtliche Daten einer Partei mit 
einer einzigen Variablen. 

Das N-Komponenten-Feld wird in Pascal mit ARRAY[1..MAX] OF 
Typ bereitgestellt, während in BASIC z.B. die Anweisung 

DIMZ (N) notwendig ist. 
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Unterprogramme werden in Pascal im Deklarationsteil als 
PROCEDURES aufgeführt (Blockstruktur), in BASIC stehen die 
SUBROUTINES hinter dem Hauptprogramm. Sie werden mit GOSUB 
Zeilennummer angesprungen und enden mit dem Rücksprungbefehl 
RETURN. 
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